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SECTION  I 


INTRODUCTION 

In  this  report  we  shall  consider  the  cruise  performance  of  a 
jet-propelled  aircraft  at  high  speed.  A characteristic  of  the  high- 
speed performance  problem  is  that  the  aerodynamic  and  engine  con- 
trols are  strongly  influenced  by  the  altitude  and  the  Mach  number. 
Hence,  in  performance  analysis,  it  is  nearly  impossible  to  obtain 
general  results  which  are  unrestrictedly  valid  for  all  flight  mis- 
sions. For  that  reason  we  shall  restrict  our  investigation  to  the 
cruise  performance  problem. 

In  general,  cruise  performance  is  characterized  by  flight  in 
the  vertical  plane  at  small  flight  path  angle  and  nearly  constant 
engine  setting.  This,  in  turn,  leads  to  slow  variations  in  the  speed 
and  altitude. 

Typical  optimum  cruise  performance  problems  involve  the  deri- 
vation of  aerodynamic  and  engine  controls  for  maximum  range  or  max- 
imum endurance  in  either  constant  altitude  flight  or  in  missions 
which  allow  variations  in  the  altitude.  The  constant  altitude  flight 
is  representative  of  commercial  flights  regulated  by  the  Federal 
Aviation  Agency.  Since  maximizing  the  range  for  a given  amount  of 
fuel  is  obviously  equivalent  to  minimizing  the  fuel  to  cover  a pre- 
scribed distance,  the  effort  of  the  author  of  this  technical  report 
is  in  line  with  the  national  policy  of  saving  energy.  On  the  other 
hand,  the  maximum  endurance  cruising  flight  is  representative  of 
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military  peacetime  missions.  Maximum  endurance  for  a given  amount 
of  fuel  is  also  equivalent  to  minimizing  the  fuel  for  the  prescribed 
peacekeeping  duration  of  a mission. 


T 


The  organization  of  this  report  is  as  follows.  After  this  in- 
troductory section,  the  equations  of  motion  and  the  problem  defin- 
ition are  discussed  in  Section  II.  Pertinent  assumptions  regarding 
the  analysis  will  be  presented  in  Section  III.  The  equations  of 
motion  using  dimensionless  variables  will  be  proposed  in  Section  IV. 
In  Section  V,  the  complete  analytical  solutions  for  the  two  prob- 
lems of  maximum  range  and  maximum  endurance  in  the  case  of  constant 
altitude  flight  will  be  presented.  The  problem  of  optimum  aero- 
dynamic and  engine  controls  for  cruise  performance  with  variations 
in  the  altitude  is  the  subject  of  the  next  two  sections.  In  Sec- 
tion VI,  an  engineering  approach  is  used  to  derive  simple  optimum 
laws  for  aerodynamic  and  engine  controls  in  both  the  maximum  range 
and  maximum  endurance  problems.  In  the  next  section.  Section  VII, 
the  same  problems  are  treated  by  optimal  control  theory  using  the 
full  set  of  equations. 
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SECTION  II 


PROBLEM  DEFINITION 


The  flight  of  a lifting,  thrusting  vehicle  in  a vertical  plane 
is  governed  by  the  equations  (Fig.  1) 


dX  „ 

— = V cos  y 


dZ  „ ^ 

— = V sin  y 


m = T cos  (e  + a) 
dt 


- D - mg  sin  y 


(1) 


mV  4^-  = T sin  (e  + a)  + L - mg  cos  y 
dt 

dm  _ c T 

dt  “ g 

The  first  two  equations  are  the  kinematic  equations.  The  next 
two  equations  are  the  force  equations  written  along  the  directions 
tangent  and  normal  to  the  flight  path.  The  last  equation  represents 
the  mass  flow  as  related  to  the  thrust.  Standard  notation  as  defined 
in  the  nomenclature  section  has  been  used.  The  flight  path  angle  y 
has  been  selected  such  that  y is  positive  when  the  velocity  is  di- 
rected above  the  local  horizontal  plane. 

The  flight  is  controlled  by  varying  the  lift  and  drag  forces, 
through  the  angle  of  attack  a,  and  the  thrust  magnitude  T.  It  is 
proposed  to  find  these  aerodynamic  and  engine  controls  to  bring  the 
vehicle  from  a certain  prescribed  initial  state  (X^  Z^, 
m^)  to  some  partially  prescribed  final  state  in  the  first  four  var- 
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iables  and  a prescribed  final  mass  rof  in  such  a way  that  a certain 
function  of  the  final  variables  is  a maximum.  In  the  problem  of 
maximum  range,  we  maximize  the  final  distance  Xf  while  in  the  prob- 
lem of  maximum  endurance  we  maximize  the  final  t^. 


W * mg 


Fig.  1.  Nomenclature. 
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SECTION  III 
ASSUMPTIONS 

To  analyze  the  system  of  equations  (1),  two  types  of  assumptions 
are  introduced.  The  first  type  consists  of  the  modeling  of  the 
medium  in  which  the  vehicle  is  flying,  namely  in  this  case  the 
earth's  atmosphere,  and  the  vehicle  aerodynamics  and  engine  charac- 
teristics. This  is  necessary  for  any  mathematical  treatment  of  a 
physical  problem.  The  second  type  of  assumptions  is  designed  to 
simplify  the  equations  of  motion  such  that  analytical  manipulation 
is  feasible.  In  this  respect  we  shall  neglect  the  terms  which  are 
small,  and  treat  as  constants  the  terms  which  are  insensitive  to 
change  along  the  flight  path.  For  the  results  to  have  any  validity, 
unnecessary  simplification  should  be  avoided  and  accuracy  of  the 
solutions  shouldn't  be  compromised. 

The  aerodynamic  characteristics  of  the  vehicle  strongly  depends 
on  the  Mach  number  defined  as 


M = — 
a 


(2) 


where  a is  the  speed  of  sound  in  the  ambient  atmosphere.  The  speed 
of  sound  is  related  to  the  atmospheric  density  p and  the  pressure  p 
through  the  relation 

a2  - JSE. 


(3) 


where  k is  the  ratio  of  the  specific  heats  of  the  ambient  air,  con- 
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sidered  as  a constant.  An  approximate  relation  between  the  pres- 
sure and  the  density  is 


where  subscript  zero  denotes  the  condition  at  sea  level  and  X and 
n are  two  constants  appropriately  selected.  Two  sets  of  values  for 
X and  n will  be  used  depending  on  the  atmospheric  layer  in  which  the 
vehicle  is  flying.  The  troposphere  is  the  first  layer  of  the  atmos- 
phere extending  from  sea  level  to  approximately  36,089  feet.  Across 
this  layer  the  atmospheric  temperature,  and  hence  the  speed  of 
sound,  varies  with  the  altitude.  Above  the  limit  of  the  troposphere, 
called  the  tropopause,  and  extending  to  approximately  66,000  feet  is 
the  stratosphere.  In  the  stratosphere  the  atmospheric  temperature 
and  the  speed  of  sound  remain  nearly  constant.  The  flight  will  be 
assumed  to  take  place  in  either  the  troposphere  or  the  stratosphere. 
For  the  sake  of  continuity  of  the  equations,  it  is  more  convenient 
to  use  the  reference  condition  at  the  tropopause.  Hence,  we  shall 
use  the  density-pressure  relationship  in  the  form 


where  subscript  * denotes  the  condition  at  the  tropopause.  The 
value  of  n is  n = 1.235  in  the  troposphere  and  n » 1 in  the  strat- 
osphere. 

The  aerodynamic  drag  and  lift  are  assumed  to  be  of  the  form 
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(6) 


D - | P S V2  CD 
L = y p S V2  CL 


where  CQ  and  are  respectively  the  drag  and  the  lift  coefficient. 

It  is  assumed  that  and  CL  are  related  by  the  parabolic  drag  polar. 


CD  +KCL 
o 


(7) 


where  the  zero-lift  drag  coefficient  and  the  induced-drag  co- 

o 

efficient  K are  functions  of  the  Mach  number 


CD  CD  (M) 

° ° (8) 

K = K (M) 

From  experimental  data,  it  is  possible  to  model  these  functions  as 
polynomials  in  M.  In  general,  it  is  assumed  that  they  are  contin- 
uous functions  in  M with  continuous  first  derivatives.  Quite  fre- 
quently, for  a given  Mach  number,  the  drag  polar,  that  is  the  plot 

of  C versus  C does  not  yield  the  minimum  of  C for  C = 0 (Fig.  2). 
L)  Li  U L 

As  shown  in  Ref.  1,  it  is  remedied  by  a simple  translation  along  the 
axis.  This  will  introduce  additional  terms  in  the  pertinent  re- 
sults. Fortunately,  in  optimum  cruise  performance,  the  operating 
point  for  the  angle  of  attack  is  near  the  point  of  maximum  lift-to- 
drag  ratio  and  not  where  the  drag  coefficient  is  minimum  (Fig.  2). 

In  this  respect,  the  accuracy  is  not  compromised  by  using  relation 

(7)  and  evaluating  CD  and  K such  that  we  have  perfect  fitting  of 

o 
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the  drag  polar  in  its  upper  part,  especially  near  the  point  of  max- 
imum lift-to-drag  ratio.  This  ratio  defined  as 


(9) 


is  a function  of  the  Mach  number  and  is  an  important  parameter  in 
performance  analysis.  If  relation  (7)  is  used,  then 


E 

max 


(10) 


It  is  assumed  that  parabolic  drag  polar  is  valid  up  to  a certain 

angle-of-attack  called  the  stalling  angle-of-attack.  At  this  point 

C reaches  a maximum  value 
L 


CL  = CL  (M)  (11) 

max  max 


and  beyond  it  the  lift  coefficient  drops  drastically  because  of  the 

turbulent  effect  due  to  flow  separation.  For  all  practical  purposes, 

we  can  consider  CL  , which  is  a function  of  the  Mach  number,  as 
max 

the  limiting  value  for  the  lift  coefficient. 

The  parameter  c is  called  the  specific  fuel  consumption.  It 
is  a function  of  the  Mach  number,  of  the  engine  throttling,  and  in 
the  neighborhood  of  the  tropopause  it  slightly  depends  on  the  al- 
titude. Generally  speaking,  for  subsonic  cruise,  there  is  a trade- 
off between  the  aerodynamic  characteristics  and  the  specific  fuel 
consumption  to  find  the  best  cruise  point.  Since  this  influence 
has  a secondary  effect  as  compared  to  the  influence  of  the  Mach 
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number  and  also  since  the  optimum  value  of  c is  usually  near  its 
minimum  value,  we  shall  assume  that  c is  minimum  and  this  minimum 
value  is  purely  a function  of  the  Mach  number. 

It  is  accurate  to  assume  that  the  acceleration  of  the  gravity 
is  constant.  We  now  turn  to  the  second  type  of  assumptions.  In 
the  equations  (1) , e is  the  angle  between  the  thrust  vector  and  the 
reference  line  from  which  we  measure  the  angle  of  attack  a.  Both 
e and  a are  small  and  the  analysis  will  be  much  simplified  by  taking 
in  the  equations 


T cos  (e  + a)  ; T 
T sin  (e  + a)  = 0 


(12) 


This  is  the  same  as  saying  that  the  thrust  vector  is  always 
aligned  with  the  velocity  vector.  The  assumption  may  not  be  valid 
in  the  case  of  flight  at  low  speed  of  a delta-wing  aircraft.  Other 
simplifying  assumptions  will  be  introduced  later  when  it  becomes 
necessary  to  use  them  in  order  to  obtain  an  explicit  solution  to  the 
problem  considered. 


SECTION  IV 


DIMENSIONLESS  EQUATIONS 


Using  the  assumptions  (12),  we  rewrite  the  equations  of  mo- 
tion using  a parabolic  drag  polar 


dX 

— = V cos  y 
at 


dZ  „ , 

— = V sin  y 


ii  ^ * T - j p S V2  (Cjj  + K C2)  - mg  sin  y 

o 

mV  = "T  P s V22  CT  - mg  cos  y 

at  L L 


dm  _ c 
dt  " g 


The  flight  is  controlled  by  the  angle-of-attack  or,  equiva- 


lently the  lift  coefficient  C^,  and  the  thiust  magnitude  T subject 
to  the  constraints 


0 ±C  < C. 

L<  Li 


0 < T < T 
— — max 


There  exists  a constraint  on  the  load  factor  n **  L/W  and  it  is  ex- 
pressed by 


n . < n < n 

min  — — max 


It  is  convenient  for  the  analysis  to  write  the  equations  of 


motion  in  dimensionless  form. 


11 


The  advantage  of  using  dimensionless  variables  is  that  if  the 
final  equations  can  be  written  in  a form  such  that  they  are  free 
of  the  physical  characteristics  of  the  vehicle,  then  the  results  of 
the  analysis  apply  to  any  aircraft  regardless  of  its  size,  mass, 
aerodynamics  and  engine  characteristics.  This  has  been  illustrated 
by  the  excellent  analysis  of  reentry  trajectories  by  Chapman  (Ref. 
2).  Chapman's  equations  are  valid  for  the  case  of  gliding  flight, 
T=0,  and  when  the  aerodynamics  characteristics  are  independent  of 
the  Mach  number.  Unfortunately,  this  is  not  true  for  our  case  and 
hence  although  we  can  put  the  equations  in  dimensionless  form  they 
still  contain  the  vehicle's  aerodynamics  and  engine  characteristics. 
Nevertheless,  the  results  obtained  are  general  enough  such  that  they 
apply  to  a large  class  of  aircraft. 

Chapman's  search  for  dimensionless  variables  has  been  by  trial 
and  error.  Here  we  shall  start  with  the  original  equations,  the 
Eqs.  (13),  and  use  a series  of  rational  and  systematic  transforma- 
tions leading  to  a logical  choice  of  the  dimensionless  variables. 

First,  since  the  aerodynamics  and  engine  characteristics  are 
functions  of  the  Mach  number,  M is  the  choice  of  the  dimensionless 
variable  to  replace  the  speed  V.  In  this  respect,  we  consider 


dV 

dt 


dM 

dt 


M dF 


(15) 


In  the  troposphere,  a is  a function  of  the  altitude,  while 
it  is  a constant  in  the  stratosphere.  Since  the  analysis  considers 


cruise  performance  in  which  the  altitude  either  remains  constant 
or  varies  slowly,  we  shall  take  da/dt  ■ 0. 

Next,  it  is  convenient  to  scale  the  weight  and  the  thrust  with 
respect  to  the  initial  weight  of  the  vehicle.  Hence,  we  can  define 
the  dimensionless  weight  w and  the  dimensionless  thrust  t as 


Now,  since  the  speed  of  sound  a will  appear  in  the  equations 
through  relation  (15),  with  relations  (3)  and  (5)  we  can  express  a 
in  terms  of  either  the  pressure  p or  the  density  p of  the  atmosphere. 
Then,  either  p or  p can  be  used  to  represent  the  altitude.  If  we 
choose  p,  then 


where  subscript  * denotes  the  condition  at  the  tropopause. 
Since 


PV‘ 


„2  2 

pM  a 


kpM 


we  can  write  the  equation  for  y 


Mw 


\a*/  8 at  2wi 


C - w cos  y 
Li 


Then  to  eliminate  kS  from  the  equation,  it  is  logical  to  define  the 
altitude  variable  as 


n 
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In  this  way,  since  p decreases  as  the  altitude  increases,  n varies 
in  the  same  direction  as  the  altitude.  The  equation  for  y is  di- 
mensionless if  the  dimensionless  time  6 is  defined  as 


With  the  dimensionless 
is 


time  so  defined. 


the  equation  for  the  range 


(a  ji  2 j 

\ /\aj 


dX 

— = M cos  y 


Consequently,  we  should  replace  the  longitudinal  distance  X and  the 
altitude  Z by  their  dimensionless  equivalents 


2 

a. 


2 

a. 


Finally,  noticing  that 


5_  /j_j  <n-1)/2n 


we  can  write  the  system  (13)  as 


dx 

de 


(l-n)/2n 


/ \ (l-n)/2n 

(t)  " 

dz  _ / n \ 
de 

*("*) 

a(t) 


cos  y 


M sin  y 


dM 

de 


dx 

de 


(n-l)/2n 


(n-l)/2n 


[t  - n_1M2(CD  + KC2)  - w sin  y] 
o 

(n  -VcL  - w cos  y) 
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The  dimensionless  specific  fuel  consumption  can  then  be  defined  as 


k 

c 


For  flight  in  the  stratosphere  (n*)  ^ = 1 and  the  equations 
are  free  of  the  parameter  n*.  To  remove  this  term  when  n j4  1,  we 
define 


B - <n„)(1-n)/2n 


(l-n)/2n 


(16) 


and  include  it  as  a coefficient  of  the  dimensionless  time  0 as  can 

be  done  in  the  equations  for  M and  y.  For  the  other  equations,  it 

requires  a new  definition  for  x,  z and  k . This  finally  leads  to 

c 

the  following  definition 


x 


, z 


Z 


w 


n 


w_ 

wi 

^i 

kpS 


L 

w. 


K = 
c 


Bg 


(17) 


In  some  analysis,  it  is  convenient  to  combine  w and  n by  using 
Miele's  dimensionless  wing  loading  defined  as  (Ref.  3). 

to  =*  * nw  (18) 

Kp  b 

With  the  dimensionless  variables  as  defined  in  Eqs.  (17)  we 


have  the  dimensionless  equations  in  their  final  form 


dx 

(l-n)/2n 

M cos  y 

de 

n 

dz 

de 

n(l-n)/2n 

M sin  y 

dM 

de  = 

ri(n-l)/2n 

w 

[t  - n-1M2( 

dl  = 
de 

r)(n-l)/2n 

wM 

(n_1«2cL  - - 

dw 

de 

— rs.i 

C 

o 

cos  y) 


First,  we  notice  that  the  two  dimensionless  variables  z and  n are 
not  independent.  Explicitly,  if  Z is  measured  from  the  tropopause 
we  use  an  exponential  of  the  form 

P = p^e~^Z  (2 

where  the  constant  6 has  different  values  in  different  atmospheric 
layers.  Then  we  have 


= eeZ 


Hence,  if  we  choose  n as  the  altitude  variable,  then  the  equation 
for  z can  be  replaced  by  the  equation  for  n 


dn  _ (n+1) / 2n 

d?  ‘ h n M sin  y 


where  the  constant  IC,  is  defined  as 
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Next,  since  Kc  is  proportional  to  c,  this  specific  fuel  consump- 
tion coefficient  is  function  of  the  Mach  number 


K = K (M) 
c c 


(24) 


It  is  assumed  that  this  function  can  be  approximated  by  a certain 
polynomial.  Finally,  the  thrust  coefficient  varies  between  its 
limits 


0 < t < T 
— — max 


(25) 


While  the  case  r = 0 is  trivial,  the  upper  limit  t 
tion  of  both  the  altitude  and  the  Mach  number. 


max 


is  func- 
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SECTION  V 


FLIGHT  AT  CONSTANT  ALTITUDE 


If  the  aircraft  is  required  to  maintain  a constant  altitude, 
Y = 0,  then  dy/d0  = 0 and  we  have  the  equilibrium  condition 


The  system  of  dimensionless  equations  (19)  is  reduced  to 


(26) 


dx  (l-n)/2n  M 
de  - n M 

dM  n(n-1)/2n 
d0  w T 


dw 

d0 

d0_ 

d0 


T 


-(n+l)/2n 
n M2 

W 


2 2 

+ ^4W-> 

M 


(27) 


where  now  n is  a constant  parameter.  The  single  control  variable 
is  the  thrust  coefficient  t subject  to  the  constraint  (25).  The 
last  of  the  Eq.  (27),  which  is  trivial,  is  introduced  to  treat  both 
the  maximum  range  and  the  maximum  endurance  problems  in  a single 
formulation.  Here,  it  should  be  emphasized  that  Kc  is  only  a func- 
tion of  the  Mach  number,  but  not  of  the  thrust,  an  assumption  which 
is  generally  valid  for  subsonic  cruise. 

Variational  Formulation 

The  problem  is  formulated  as  an  optimal  control  problem.  The 
equations  (27)  are  the  state  equations.  The  initial  conditions  are 
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e - 0,  X - 0,  M = Mi,  w = w1  = 1 (28-a) 


It  is  required  to  find  the  thrust  magnitude  control  x,  and  sub- 
sequently the  angle-of-attack  to  maintain  level  flight,  such  that 
at  the  final  time,  for  a prescribed  and 

0 - 6f,  x - xf,  M - Mf,  w - wf  (28-b) 

For  the  maximum  range  problem,  we  maximize  x^  while  allowing 
the  final  time  0^  to  be  free,  and  for  the  maximum  endurance  problem, 
the  final  time  is  to  be  maximized  while  the  range  xf  is  free. 

Using  the  maximum  principle,  we  introduce  the  adjoint  vector 
p = (p  ,p„,p  ,pQ)  to  form  the  Hamiltonian 

X PI  W b 

/i  \/o  -(n+l)/2n  „ 2 2 

H = Pe  + n(  “ }/  n Mpx  - U - M PM(CD  + Mr_)  + 

o M 

t , (n-l)/2n  „ . 

— (n  pu  - wK  p ) 

w rM  c w 


(29) 


where  the  components  of  p are  governed  by  the  equations 


dp 

d0 

dp 

d0 


M 


3H  (l-n)/2n  n (n+1)/2n  2 . 

3M  " " n Px  + w PM  (2MCD  + M CD  + 

o o 


k'p2w2  2Kn2w2  > 

M2  ' M3  <=  " 


(30) 


£ - - ft  - (A  + ) + P„, 

' w M • w 
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de 


3H 

36 


- 0 


where  the  prime  denotes  derivatives  taken  with  respect  to  M. 

The  problem  is  solved  by  integrating  the  state  and  adjoint 
equations,  subject  to  the  end  conditions  (28-a)  and  (28-b)  while 
selecting,  at  each  instant  6,  the  thrust  control  t such  that  the 
Hamiltonian  H is  an  absolute  maximum. 

Optimum  Aerodynamics  and  Engine  Thrust  Controls 
The  problem  has  a number  of  integrals 


(31) 


where  and  are  constants  of  integration. 

In  the  maximum  range  problem,  if  the  final  time  is  free. 


In  the  maximum  endurance  problem,  if  the  final  distance  is 
free,  = 0. 

To  select  the  optimum  thrust  control,  we  consider  the  function 

e = n(n"1)/2npM  - Kcwpw  (32) 

called  the  switching  function.  Then,  to  maximize  H, 

If  6 > 0 we  choose  t = t 

max 

If  ® < 0 we  choose  t ■ 0 
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If  © * 0 for  a finite  time  interval  6e [e^^]  we  choose  r = var- 
iable. The  optimum  trajectory  is  a combination  of  subarcs: 

Maximum  thrust  arc  in  which  t ■ t 

max 

Coasting  arc  in  which  t = 0 
Sustaining  arc  in  which  t = variable 


It  is  assumed  that  maximum  thrust  arc  and  coasting  arc  for  both  the 
maximum  range  and  maximum  endurance  problems  are  of  short  duration. 
This  rather  intuitive  assessment  will  be  confirmed  later  by  rigor- 
ous analysis. 

Along  the  main  portion  of  sustaining  flight,  we  have  for  all 
time  0 


(n-l)/2n 
n PM 


K wp 
c w 


(33) 


By  taking  the  derivative  of  this  equation  with  respect  to  0, 


(n-l)/2n  dpM 


d0 


v dw,„^  dM  Pw 

Kp  -J7T+K  wp  "Tr  + K W -rr— 
crw  d0  c rw  d0  c d0 


(34) 


Using  the  state  and  adjoint  equations,  together  with  relation 
(33)  in  Eq.  (34),  we  have  the  condition  for  optimum  sustaining 


flight 


- M2PM 


(Cn  [2  + n(1  n>/2nMK  + (CD  kc)h1 
o o 


(35) 


t2  + n(1  n,/2nMK  - (KK  )„]} 
m4  c c M 


where,  following  Miele  (Ref.  3),  to  simplify  the  notation,  we  have 
defined  the  logarithmic  derivatives 
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d(Log  CD  ) 

o 

d(Log  M) 


d (Log  K)  _ M dK 
N d (Log  M)  K dM 


d(Log  Kc)  M dKc 
CCM  ' d(Log  M)  = Kc  dM 


On  the  other  hand,  we  can  write  the  Hamiltonian  integral  H = 0,  for 
the  case  of  sustaining  flight  as 

2 

+ ( +Kn^)  (37 

i / w u „** 

o M 

By  eliminating  the  adjoint  between  the  Eqs.  (35)  and  (37) 

we  have  the  relation  between  the  weight  w and  the  Mach  number  M for 

the  general  case  of  sustaining  flight  at  constant  altitude. 

For  the  case  of  maximum  range,  free  final  time,  we  have  = 

0.  The  elimination  of  pu  leads  to  relation 

M 


U)  * M 


'l  + n(1'n)/2nMKc  + (CD  Kc)m 

O 

3 + - (KKc)m 


For  the  case  of  maximum  endurance,  free  range,  we  have  C ^ 


This  gives 


A + na'n)/2V  + ccn  k ) 


GO  * M 


D *cM 


K 2 + na^)/Z"MK  - CKKc)m 
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For  a prescribed  flight  altitude,  Eq.  (38),  upon  solving,  gives  the 
optimum  Mach  number  as  function  of  the  weight  in  the  case  of  maxi- 
mum range.  Likewise,  Eq.  (39)  gives  the  Mach  number  as  function  of 
the  weight  in  the  case  of  maximum  endurance. 

To  assess  the  influence  of  the  altitude  and  weight  on  the  Mach 
number,  we  define  the  pressure  ratio 


n n* 

6 = -2- 

P*  n 


(40) 


and  the  specific  fuel  consumption  coefficient 


k n*(1"n)/2nK 

eg*  c 


(41) 


The  Eqs.  (38)  and  (39)  have  their  final  form: 
For  maximum  range 


a)  = M 


Cn  , /I  + 6(n_1)/2nMk_  + (C„  kj 
2 1/  *o 


D c M 
o 


3 + 6(n-1)/2nMk  - (Kk  )M 

C C M 


(42) 


For  maximum  endurance 


a)  = M 


'l  + 6(tl_1)/2nMk  + (C_  k )M 

c D c M 
o 

2 + <5(n-1)/2nMk  - (Kk  )M 

C CM 


(43) 


If  point  performance  analysis  is  considered,  by  maximizing  at 
each  instant  the  rate  of  change  of  the  distance  for  the  maximum 
range  problem,  and  minimizing  the  mass  flow  for  the  endurance  prob- 
lem, we  shall  come  up  with  two  equations  identical  to  the  Eqs,  (42) 
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and  (43)  without  the  term  This  term  is  the  correc- 

c 

tional  term  due  to  the  slow  decrease  of  the  Mach  number.  In  gen- 
eral, it  is  a small  term. 

The  lift  coefficient  is  obtained  from  the  equilibrium  con- 
dition, Eq.  (26),  written  in  terms  of  the  Mach  number  as 

ui  « M2Cl  (44) 


Hence,  from  Eqs.  (42)  and  (43),  we  have: 
For  maximum  range 


1 + 6(n~1)/2nMk  + (C  K ) 

c , c*  V/  ! “olH 

L L 3 + 6(n“1)/2nMk  - (Kk  ) 

c c M 


(45) 


For  maximum  endurance 


/2  + + (CD<>kc)M 

L L 2 + 6<n-1)/2nMk  + (Kk  ) 

C CM 


(46) 


where 


* 


(47) 


is  the  lift  coefficient  corresponding  to  maximum  lift-to-drag  ratio. 

In  terms  of  the  lift-to-drag  ratio,  we  have: 

For  maximum  range 

(48) 


C 

C 


L 

D 


E 

max 


[l+6(n_1)2nMk  +(C  k )M]1/2[3+6(n''1)/2nMk  -(Kk  )M]1/2 

C U C M C C M 

O 

2+6 (n'1) /2nMkc+(l/2)(CD  /K)M 

o 
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For  maximum  endurance 


[2+6(n"1)/2nMk  +(C  k ) ]1/2[2+«(n  1)/2nMk  -(Kk  ) ]1/2 

C C M C C M 

CD  maX  2+6(n_1)/2nMk  +(1/2)(C  /K)m 

C DM 

O 

From  an  operational  point  of  view,  if  the  instantaneous  Mach 
number  is  used  as  control  variable,  then  it  is  obtained  from  either 
Eq.  (42)  or  Eq.  (43),  depending  on  the  case.  The  thrust  is  then 
adjusted  to  maintain  constant  altitude  flight.  Analytically,  since 
the  Mach  number  varies  slowly  along  the  flight,  approximate  value 
for  the  thrust  is  T=D. 

This  gives 


nr  = M2(C  + *f-) 
o M 


Hence,  for  the  case  of  maximum  range 

2+6(n-1)/2nMk  +(1/2) (n  /K)u 
2T_  m C Do  M *2 

kpS  3+6(n-1)/2nMk  -(Kk  ) D 

C CM 


For  the  case  of  maximum  endurance,  we  have  the  expression  for 


the  thrust 


2+6<"-1)/2nMk+(l/2»(C  /K) 

2T_  . C °o  " *2 

kps  (n-D/an^  D 

c c M 


where 


CD  “2CD 
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is  the  drag  coefficient  corresponding  to  maximum  lift-to-drag  ratio. 

If  exact  variational  analysis  is  considered,  then  along  the 
optimum  trajectory  for  sustaining  flight,  the  Eq.  (42),  or  (43)  must 
be  constantly  satisfied,  depending  on  the  case. 

Hence,  we  write 


with 


l+6(n-1)/2nMkc+(CDokc)M 

^ 3+6(n_1)/2nMk  -(Kk  ) 

c c M 


(54) 


(55) 


for  the  case  of  maximum  range  and 


2+6(n-1)/2nMkc+(CD^kc)M 
E 2+<5(n-l)/2n>flt  _(Kk  , 

C C M 


(56) 


for  the  case  of  maximum  endurance. 

By  taking  the  derivative  of  Eq.  (54)  with  respect  to  0,  the 
variable  thrust  control  will  appear  linearly  and,  upon  solving,  we 
have 


2T 

kpS 


4+(Acd  /K)m 

o 

4+ (AC  /K)u+26  (n-l)/2v 

DM  C 


(l+A) 

2 


* 2 
cdm 


(57) 


The  bracketed  coefficient  Is  the  correctional  factor  to  either 
Eq.  (51)  or  (52).  The  effect  is  due  to  the  term  26  (n_1^2nMkc. 
Because  of  this  term,  the  thrust  is  slightly  less  than  the  drag, 
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hence  accounting  for  a slow  decrease  of  the  Mach  number  along  the 
flight  path  in  both  the  maximum  endurance  and  maximum  range  cases. 
Subsonic  and  Supersonic  Flight  Domains 

Since  n is  constant,  it  is  convenient  to  represent  the  flight 
domain  in  the  (M,w)  space.  In  this  space,  we  plot  the  relation 
(54)  which  represents  the  trajectory  along  the  sustaining  flight 
arc,  with  A = for  the  case  of  maximum  range  and  A = A£  for  the 
case  of  maximum  endurance. 

In  this  report,  two  typical  military  aircraft  are  used  for  the 

numerical  application.  They  are  the  transport-tanker  KC-135A  and 

the  fighter  F-4C.  Their  aerodynamic  characteristics  are  given  in 

Appendix  A.  For  both  aircrafts,  the  value  k is  very  small  and  is 

c 

practically  independent  of  the  Mach  number.  Hence,  in  the  (M,<d) 
space  we  have  the  simplified  equations  for  the  trajectory  in  sus- 
tained flight 


(o  - M 


1 + C_ 


JM 


3 - 


(58) 


for  the  case  of  maximum  range,  and 

..sJsJTs 

“ k 2 - Km 


(59) 


for  the  case  of  maximum  endurance. 

In  Appendix  A,  for  each  aircraft,  the  functions  Cp  and  K are 

o 

given  in  terms  of  the  Mach  number.  The  Eqs.  (58)  and  (59)  are 
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plotted  in  Fig.  3 for  the  KC-135A  and  in  Fig.  4 for  the  F-4C. 


At  low  Mach  number,  the  influence  of  the  Mach  number  on  aero- 
dynamic characteristics  is  negligible  and  the  curves  are  nearly 
the  parabolas 


<i) 


* 

— M^ 
i/3 


* 2 

u)  = CL  M 


(60) 


where  C is  the  asymptotic  constant  value  of  C for  low  Mach  num- 
L L 

ber.  Hence,  at  low  Mach  number,  maximum  endurance  flight  is  con- 

* 

ducted  at  * CL  while  maximum  range  flight  is  conducted  at  = 
Cl//3.  At  high  Mach  number,  as  the  logarithmic  derivative 
approaches  the  value  2 for  the  case  of  maximum  endurance  and  the 
value  3 for  the  case  of  maximum  range,  the  curves  tend  respectively 
to  their  vertical  asymptotes.  From  the  figures,  it  is  seen  that  in 
both  cases,  optimum  Mach  number  has  an  upper  limit  and  this  upper 
limit  is  in  the  subsonic  range.  Also  for  the  same  w = 2W/kpS,  the 

Mach  number  for  maximum  range  is  higher  than  the  Mach  number  for 

* 

maximum  endurance.  At  low  speed,  since  is  the  same,  it  is  seen 
from  Eqs.  (60)  that  their  ratio  is  M^/M^.  = (3)^^  = 1.316.  This 
ratio  is  smaller  for  higher  id  which  corresponds  to  the  case  of  sus- 
tained flight  at  high  altitude  and/or  high  wing  loading. 

The  curves  can  be  used  if  they  are  within  the  flight  domain 
of  the  aircraft.  More  specifically,  since  the  lift  coefficient 

is  bounded  by  the  maximum  lift  coefficient  (M),  which  is  a 

max 
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function  of  the  Mach  number,  in  the  (M,io)  space,  the  domain  of 
flight,  also  called  the  flight  envelope,  is  bounded  by  the  curve 

(o  = CL  (M)  M2  (61) 

max 

At  low  Mach  number,  is  nearly  constant,  and  the  curve  is 

max 

approximately  a parabola.  At  high  Mach  number,  depends 

max 

strongly  on  the  Mach  number,  and  we  have  two  cases: 

For  subsonic  aircraft  configuration,  when  the  Mach  number  is 

approaching  unity,  CL  drops  drastically  and  we  have  the  case  of 

max 

Fig.  3.  In  practice,  the  domain  of  flight  is  further  limited  by 
the  buffeting  effect.  At  high  angle-of-attack,  before  the  lift  co- 
efficient reaches  its  maximum  value  C , appreciable  flow  separ- 

max 

ation  from  the  wing  occurs  at  a certain  lift  coefficient  CT  < CL 

B max 

called  the  buffeting  coefficient.  Beyond  this  lift  coefficient, 

the  aircraft  begins  to  shake  due  to  the  turbulence  associated  with 

flow  separation.  Like  C , the  buffeting  lift  coefficient  C 

L L_ 

max  B 

is  a function  of  the  Mach  number.  Hence,  if  we  plot  Eq.  (61)  with 

C.  replaced  by  C , we  have  another  curve  delimiting  the  buffet- 
L L_ 

max  B 

ing  region. 

In  any  instant  of  sustained  flight,  for  any  prescribed  w, 
there  are  4 Mach  numbers  of  interest  namely  M^,  M2 
3).  The  lower  bound  Mach  number  M^  is  the  stalling  Mach  number  while 
the  upper  bound  Mach  number  M2  is  the  critical  Mach  number  due  to 
the  subsonic  nature  of  the  aircraft  aerodynamic  configuration.  M^ 
and  M2  are  coincident  when  the  function  (61)  has  its  maximum  value. 


. Mj;  and  Mr  (Fig. 
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The  corresponding  Mach  number  is  obtained  by  writing  the  condition 

2 

for  the  maximum  value  of  the  function  M which,  in  terms  of 

max 

the  logarithmic  derivative,  can  be  written  as 


<CT 


max 


) +2  = 0 
M 


(62) 


where,  as  usual,  subscript  M denotes  the  logarithmic  derivative. 

The  resulting  value  of  w corresponds  to  the  aerodynamic  ceiling. 

The  Mach  number  M^  is  the  Mach  number  for  maximum  range  and 
the  Mach  number  M^,  is  the  Mach  number  for  maximum  endurance.  They 
are  obtained  by  solving  the  Eq.  (54)  with  A = A^  and  A = A^  respec- 
tively. 

For  supersonic  aircraft  configuration,  the  function  (M) 

max 

decreases  when  the  Mach  number  approaches  unity  and  then  tends  to 
a near  constant  value,  lower  than  the  constant  value  for  low  sub- 
sonic regime.  Hence,  in  the  supersonic  regime,  the  function  (61) 
is  again  near  parabolic  as  seen  in  Fig.  4.  In  this  case,  while 
the  lower  bound  Mach  number  M^  is  still  the  stalling  Mach  number, 
the  upper  bound  Mach  number  M2  is  either  the  critical  Mach  number 
due  to  the  effect  of  compressible  fluid  or  the  maximum  Mach  number 
that  maximum  thrust  can  deliver  in  level  flight.  This  Mach  number 
is  obtained  by  solving  the  equation 


(63) 
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Fig.  3.  Domain  of  Flight  and  Optimum  Trajectories  for  Subsonic  Aircraft. 


Fig.  4.  Domain  of  Flight  and  Optimum  Trajectories  for  Supersonic  Aircraft. 


Optimum  Switching 

The  general  problem  of  constant-altitude  flight  is  formulated 
as  follows.  At  a given  constant  altitude,  the  initial  and  the  final 
Mach  numbers  M^  and  are  prescribed.  For  a given  initial  weight 
and  a final  weight  W^,  find  the  optimum  aerodynamics  and  engine 
controls  to  maximize  the  range  or  the  endurance. 

To  solve  the  problem,  we  use  the  dimensionless  wing  loading 
u>  to  replace  the  weight.  For  any  given  aircraft  aerodynamic  config- 
uration and  engine  characteristic  we  plot  the  optimum  relation  (54) 
between  u>  and  M in  the  (M,u>)  space.  (Fig.  5).  This  curve  repre- 
sents the  optimum  subarc  for  variable  thrust  sustaining  flight. 

In  general,  a complete  optimum  trajectory  connecting  the  initial 
point  (M^.oj^)  and  the  final  point  (M^.to^),  inside  the  permissible 
domain  of  flight,  includes  an  initial  arc  which  can  be  a null 
thrust  arc,  x = 0,  or  a maximum  thrust  arc  x = x joining  the  in- 
itial  point  and  the  variable  thrust  arc,  and  also  a final  arc  with 
x = 0,  or  x = x leading  to  the  final  point. 

The  figure  represents  the  case  where  the  trajectory  starts 
with  a maximum  thrust  arc  to  increase  the  Mach  number  to  a value 
such  that  the  optimum  relation  (54)  is  satisfied.  From  that  point, 
the  optimum  Mach  number  is  maintained  such  that  the  relation  is 
identically  satisfied.  If  the  final  Mach  number  is  free  or  at  a 
value  to  the  left  of  the  curve,  variable  thrust  is  maintained  until 
w reaches  At  that  point,  the  thrust  is  reduced  to  zero  and  the 

aircraft  coast  flies  to  the  final  value  M^  or  to  the  stalling  Mach 
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number,  providing  either  the  maximum  range,  or  the  maximum  endur- 
ance depending  on  the  expression  used  for  the  function  A.  Along 
the  sustaining  arc  the  optimum  lift  coefficient  is  given  by  either 
Eq.  (45)  or  (46),  and  the  optimum  modulated  thrust  is  given  by  Eq. 
(57).  Along  the  initial  maximum  thrust  arc,  x = x^^  and  along  the 
final  coasting  arc  x = 0,  while  the  lift  coefficient  is  adjusted 
such  that  the  equilibrium  condition  (26)  is  satisfied.  It  appears 
now  that  the  maximum  thrust  arc  and  the  null  thrust  arc  are  short 
arcs  necessary  to  connect  the  boundary  Mach  number  to  or  from  the 
optimum  Mach  number  along  the  sustaining  arc. 

The  junction  of  two  arcs  of  different  types  is  called  a 
switching.  The  direction  of  optimum  switching  can  be  decided  upon 
by  the  following  rule  (Ref.  4).  Let 


H 

H 


1 

2 


(H>x 

0i)T 


o 


x 

max 


and  define  the  switching  function  0 as 


(64) 


0 = H2  - Hx  (65) 

Then,  the  switching  is  from  x * 0 to  x = x^^  if  and  only  if  at  the 
switching  point 


d_ 

d0 


0 


> 0 


(66) 


where  in  taking  the  derivatives,  the  state  and  adjoint  equations 
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Fig.  5.  Optimum  Switchings. 


are  used  with  t = 0.  If  the  inequality  reverses,  then  the  switch- 
ing is  from  t = t to  t = 0. 

° max 

In  applying  the  theorem,  using  the  expression  (29)  for  the 
Hamiltonian,  we  have 


max  - (n-l)/2n  „ . 

(n  pu  - K wp 

w Me  w 


(67) 


By  taking  the  derivative  of  this  equation  with  respect  to 

the  dimensionless  time  6,  keeping  in  mind  that  x is  function  of 

max 

M,  hence  is  time  varying,  but  also  noticing  that  the  condition  (66) 
is  evaluated  at  the  time  of  the  switching,  hence  8 = 0,  we  have 
the  condition  for  switching  from  x = 0 to  x = x 

max 


Mp, 


M 


{CD  [2  + 6(n-1)/2nMkc  + (CD  kc)M]  - *f-  12  + 6(n"1)/2nMkc 
o o M 


" (Kkc)M]} 


c2n 


l/n 


(68) 


For  the  case  of  maximum  endurance,  = 0,  > 0. 

switching  the  Hamiltonian  integral  is  reduced  to 


-(n+l)/2n  „ ,2 

M VS  + " C1  > ° 

o M 


Hence,  pw  > 0 and  the  condition  (68)  leads  to 
M 


w < M 


Do  1/2 

IT 


At  the  time  of 


(69) 


(70) 


Switching  from  x = 0 to  t = t is  conducted  in  the  region  to  the 

max 

right  of  the  sustaining  trajectory  and  from  x = xmax  to  x * 0 in 
the  region  to  the  left  of  the  curve. 
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For  the  case  of  maximum  range,  = 0,  > 0. 

The  Hamiltonian  integral  is  reduced  to 


MPM  (r  , KjA  _ 1/n  . n 

— (CD  + — ) - C2n  >0 

o M 


(71) 


Hence,  we  also  have  PM  > 0,  and  combining  the  Eqs.  (68)  and  (71) 


we  have  the  condition  for  a switching  from  x = 0 to  x = x 


max 


a)  < M 


D°  Al/2 

r ar 


(72) 


We  have  the  same  conclusions  for  the  different  regions  of  optimum 
switching  (Fig.  5). 

The  sustaining  arc,  as  given  by  Eq.  (54)  divides  the  admissible 
flight  domain  in  the  (M,w)  space  into  two  regions  I and  II.  In  the 
region  I,  the  switching  is  from  x = xmax  to  x = 0,  while  in  the  re- 
gion II  the  switching  is  from  x = 0 to  x = Tmax*  Let 

N:  coast  arc,  null  thrust,  x *=  0 

M:  maximum  thrust  arc,  x = x 

max 

S:  sustaining  arc,  x = variable 


If  the  initial  point  is  in  region  I,  the  trajectory  starts 
with  maximum  thrust  arc,  while  it  starts  with  a coast  arc  if  the 
initial  point  is  in  region  II.  Likewise,  if  the  final  point  is  in 
the  region  I,  the  trajectory  terminates  with  a coast  arc  while  it 
terminates  with  a maximum  thrust  arc  if  the  final  point  is  in  the 
region  II.  If  the  final  Mach  number  is  free,  at  the  final  time 
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Pw  = 0,  and  by  the  definition  (32),  © is  necessarily  negative. 

Mf 

Hence,  the  optimum  trajectory  terminates  with  a coasting  arc. 

With  these  considerations  we  have  six  types  of  optimum  trajec- 
tories as  shown  in  Fig.  6.  In  general,  for  both  the  maximum  range 
problem  and  the  maximum  endurance  problem,  the  weight  of  the  fuel 
W^  - W^  is  large  enough  such  that  the  main  portion  of  the  optimum 

i 

trajectory  is  a sustaining  arc.  We  have  ruled  out  the  case  of  pure 
coasting  flight  AW  = 0,  and  the  case  of  pure  maximum  thrust  flight 
since  they  are  special  cases  with  the  sustaining  arc  vanishing. 
Using  the  theory  in  Ref.  4,  it  can  be  proved  rigorously  that  there 
are  no  more  than  two  switchings  and  Fig.  6 summarizes  all  possible 
optimum  trajectories. 

j 
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SECTION  VI 


CRUISE  PERFORMANCE 


Again,  we  consider  the  problem  of  maximum  range  and  the  prob- 
lem of  maximum  endurance.  While  in  the  previous  section  the  flight 
is  constrained  to  perform  at  constant  altitude,  this  restriction 
is  removed  in  this  section.  Obviously,  performance  is  improved 
without  constraint  on  phase  space.  This  is  one  of  the  most  attrac- 
tive aspects  of  cruise  performance. 

Generally  speaking,  optimum  cruise  is  performed  at  very  small 
flight  path  angle.  Hence,  we  can  take  approximately  cos  y ~ 1 and 
sin  y ' y and  write  the  general  equations  in  the  form 

dx  _ (1-n) /2rL, 

de  n ™ 


dn  _ K (n+i)/2n^ 

de  ‘ V 


dM  n(n~1)/2n 
d8  w 

(n-l)/2n 
dy  = n 

d0  wM 


[T  - n VcCp  + kc2)] 

o 

(n_1M2cL  - W) 


(73) 


dw 

de 

de 

de 


- K T 
c 

1 


where,  in  the  equation  for  M,  we  have  neglected  the  small  component 
of  the  gravity  w sin  y as  compared  to  the  drag  force. 

The  optimum  aerodynamics  and  engine  controls  for  both  problems 
will  be  derived  by  two  different,  but  closely  related,  analytical 
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approaches.  The  first  one  will  be  given  in  this  section. 


r r 


Equilibrium  Cruise 

For  cruise  performance,  the  flight  path  angle  not  only  remains 
small  but  varies  slowly  along  the  cruise  portion  of  the  trajectory. 
Hence,  we  can  take  dy/d0  = 0,  and  obtain  the  equilibrium  cruise 
condition 


(74) 


This  has  the  result  of  replacing  the  lift  coefficient  by 
the  altitude  n as  one  of  the  controls,  the  other  being  the  dimen- 
sionless thrust  t. 

On  the  other  hand,  since  the  altitude  varies  slowly,  we  can 
remove  the  equation  for  n from  consideration.  Both  assumptions  are 
justified  by  the  fact  that  for  most  aircraft  in  use,  over  a range 
of  several  hundred  nautical  miles  the  gain  in  altitude  in  optimum 
cruise  is  only  a few  thousand  feet.  Hence,  both  the  variations  in 
n and  in  y are  negligibly  small. 

With  these  assumptions,  the  system  (73)  is  reduced  to 


dx 

d0 

dM 

d0 

dw 

d0 

d0 

d0 


n(1"n)2nM 


n(n-l)/2n  -12  KqV 

It  - n M (CD  + )] 

o M 


- K T 
c 


(75) 


41 


which  is  identical  to  system  (27)  for  level  flight.  Nevertheless, 
the  main  difference  here  is  that  n is  not  a constant  parameter,  but 
in  this  formulation,  it  is  considered  a control  variable  which  is 
varying  slowly  such  that  approximately  dn/d0  = 0.  It  should  be 
noted  that  the  assumptions  are  introduced  solely  to  obtain  an  anal- 
ytical solution  for  the  optimum  aerodynamics  and  engine  controls. 
Once  the  solution  for  the  controls  has  been  obtained,  the  elements 


of  the  trajectory,  and  in  particular  the  altitude  q and  the  flight 
path  angle  y are  obtained  by  integrating  the  full  set  of  equations, 
the  Eqs.  (73),  or  better,  the  exact  system  (19)  with  the  equation 
in  z replaced  by  the  equation  (22)  in  q. 

The  analysis  is  the  same  as  in  the  case  of  level  flight.  In 
particular,  along  the  main  portion  of  sustaining  flight  we  have  the 
relations  (35)  and  (37)  rewritten  here  for  convenience 


= {CD  [2  + n(1_n)2nMKc  + (CD  Kc)m]  - 

o o 


Kf|/-VL  [2  + n(1-n)/2nMK  - (KK  )M]} 

w4  c c M 


2 

_ (n+l)/2n  1/rr.  _ M PM  KqV 

Cln  + °2n  * — (CD  + m4  > 

O M 


Again,  we  consider  the  expression  for  the  Hamiltonian 


H “ C 


(1-n)  /2ir.  n-(n+l)/2n  - KqV 

1 + C2n  >1 M pM(CD  + — — ) + 


o M 


i !"<n'1)/2\,  - Kc”pw) 
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With  respect  to  n»  considered  as  a control,  H is  maximized  when 
dH/dn  * 0. 

We  have 

0 

2 2 

i_ 

Since  the  Mach  number  varies  slowly  the  thrust  is  nearly  equal  to 
the  drag  and  we  have 


M2p 


(l-n)C2n1/nM  + ~ [(n+l)CD  - (3n-l)^H  + (n-l)^  = 0 


tit  = M2(Cd  + KC2) 
o 

Substituting  into  the  equation  above,  we  have  the  relation 

. . 2nM2pu  „ 2 2 

(n-l)C2n1/nM  = — (CD  - (78) 

o M 

2 

By  eliminating  M pM/w  among  the  Eqs.  (76),  (77)  and  (78)  we  have 
two  relationships  for  the  weight  w,  the  altitude  n and  the  Mach  num- 
ber M along  the  sustaining  arc. 

The  following  particular  cases  are  of  interest: 

Maximum  Range 

In  the  problem  of  maximum  range  with  free  final  time,  we  have 
Cx  = 0.  Then  the  Eqs.  (76),  (77)  and  (78)  become 

mV 


C2«17"m  ' T®  (CD  12  + nU‘n>/2"MK_  + <C„  K.)„]  - 
o 


c ' D c M 
o 


~ [2  + n(1_n)/2nMKc  - (KKc)Mn 
M 


c//'V -^<cD  + 


43 


. , 2nM2p„  „ 2 

<-1>c!"  ""  ■ <cD  -7-) 

o M 


Eliminating  C ^ between  the  last  two  equations,  we  have 

n + 1 2 


(79) 


On  the  other  hand,  the  elimination  of  between  the  first  two  equa- 
tions yields 


2 h 

-M  2v-^ 


where 


If  4“ 


! + 6(n"1)/2nMk  + (C  k ) 

c D c M 

o 

3 + 6(n_1)/2nMk  - (Kk  )M 
c c M 


(80) 


(81) 


where  6 = p/p*  is  the  pressure  ratio.  The  equilibrium  condition 
(74)  and  the  optimum  relations  (79)  and  (80)  provide  the  optimum 
lift  control  law,  the  optimum  Mach  number  and  the  optimum  cruise 
altitude  for  the  problem  of  maximum  range. 

Combining  the  Eqs.  (74)  and  (79),  we  have  the  optimum  lift  co- 
efficient as  given  in  Ref.  1 


V 3n  - 1 


(82) 


In  terms  of  the  lift-to-drag  ratio,  we  have 


✓(n+1) (3n-l) 
2n 


max 


(83) 
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Thus,  in  the  stratosphere  where  n = 1,  the  flight  is  at  maximum 
lift-to-drag  ratio,  while  in  the  troposphere  where  n = 1.235 


= 0.99546  E 
Cp  max 


(84) 


In  terms  of  the  aerodynamic  control,  optimum  cruise  for  max- 
imum range  is  quite  different  than  in  the  case  of  level  flight. 

This  can  be  interpreted  as  while  holding  down  the  aircraft  to  main- 
tain constant  altitude,  performance  is  diminished,  but  if  the  re- 
striction is  removed,  the  aircraft  will  climb  to  the  optimum  al- 
titude and  cruise  in  the  optimum  balance  between  lift  and  drag. 

Eliminating  w between  the  Eqs . (79)  and  (80),  we  have  the  equa- 
tion for  evaluating  the  optimum  Mach  number 


Ar  = (n+1) / (3n-l) 


(85) 


Explicitly,  we  have 


2(n-l)6(n  1)/2nMkc+(n-l)(Kkc)M+(3n-l)(CD  k.)M  = 4 (86) 

o 

If  the  flight  is  performed  in  the  stratosphere,  n=l,  and  we 


have 


= 1 


(87) 


\ ma  / M 


It  is  seen  that,  the  optimum  Mach  number  is  function  of  the 
specific  fuel  consumption  coefficient  k^  defined  by  Eq.  (41)  and 
the  maximum  lift-to-drag  ratio.  Hence,  while  for  flight  in  the 
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troposphere,  the  optimum  Mach  number,  as  given  by  Eq.  (86),  slightly 
depends  on  the  altitude  through  the  pressure  ratio  <5  besides  its 
aerodynamics  and  engine  characteristics  dependence,  for  flight  in 


the  stratosphere,  to  the  accuracy  of  the  assumptions  introduced 
above,  the  optimum  Mach  number  as  given  by  Eq.  (87)  is  a constant. 

For  flight  in  the  stratosphere  we  have  the  following  pertinent 
remarks : 

1.  The  function 


k 

g-5-  - f (M)  (88) 

max 

is  the  similarity  function  for  cruise  performance.  That  is  to  say 
if  two  aircrafts  have  the  same  function  f(M),  they  will  fly  at  the 
same  Mach  number  for  maximum  range, 

2.  Since  the  Mach  number  is  constant,  and  K are  constant 

o 

and  by  Eq.  (79)  ui  = nw  is  constant.  Since  w decreases,  n increases 
and  hence  the  altitude  increases. 

3.  As  mentioned  above,  the  optimum  law  is  only  approximate. 
Using  this  approximate  law  in  the  exact  equations  for  integration 
the  resulting  product  qw  is  not  a constant  but  is  slowly  varying. 

4.  With  the  approximation  u>  = constant,  we  write 


(89) 


Hence,  if  an  exponential  atmosphere  is  used,  the  gain  in  al- 
titude over  the  whole  sustaining  arc  is 
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and  is  of  the  order  of  a few  thousand  feet.  Over  a range  of  a few 
hundred  nautical  miles  and  considering  the  trajectory  as  nearly  a 
straight  line,  the  flight  path  angle,  y,  expressed  in  radians  is 
negligibly  small. 

A complete  optimum  trajectory  will  include  a climb  to  a cer- 
tain optimum  altitude  where  the  variable  thrust  phase  will  begin. 

As  in  the  case  of  level  flight,  the  optimum  trajectory  ends  with 
a climb  to  the  final  prescribed  altitude  with  the  final  prescribed 
Mach  number,  or  a glide  to  the  final  altitude.  This  complete  prob- 
lem is  difficult  to  solve  analytically. 

Approximately,  we  can  consider  a suboptimum  climb  to  an  altitude 
near  the  tropopause  with  a resulting  weight  to  be  used  as  the  in- 
itial weight  for  the  sustaining  phase.  Then  the  optimum  cruise  al- 
titude and  the  optimum  Mach  number  are  obtained  by  solving  the  two 
equations  (79)  and  (86)  for  flight  in  the  troposphere.  For  flight 
in  the  stratosphere,  the  Mach  number  is  obtained  from  Eq.  (87)  and 
then  the  altitude  is  given  directly  by  Eq.  (79)  with  oj  = nw^. 

The  equation  for  the  optimum  Mach  number  can  be  solved  only 
for  given  aerodynamics  and  engine  characteristics  as  will  be  car- 
ried out  later  in  this  section.  Nevertheless,  following  Miller 
and  Koch  (Ref.  1),  we  can  generate  some  useful  graphs  for  use  with 
any  aircraft.  In  this  purpose,  we  rewrite  Eq.  (79)  as 
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(91) 


or 


W _ M2 

kpSCDo  2/(3n-l)KCD 


I (— 

V n+1  V kp . 

K L/ 

O 

where  6 = p/p*  is  the  pressure  ratio, 
of  the  two  functions 


W V 

SC-.  )&  = 


E M 
max 


(92) 

Then,  we  consider  the  graphs 


and 


y = 


E 

max 


y 


ft 

6 


where  the  parameter  ft  is  defined  as 


(93) 


(94) 


n = 


(95) 


o 

Fig.  7 presents  the  plot  of  Eq.  (93)  for  different  values  of  E , 

max 

and  the  plot  of  Eq.  (94)  for  different  values  of  the  modified  wing 
loading  ft.  For  convenience,  in  this  plot,  the  pressure  ratio  6 has 
been  expressed  in  terms  of  the  altitude,  and  the  ordinate  is  plotted 
in  terms  of  Log  y. 

As  application,  consider  a flight  in  the  stratosphere.  For  a 
certain  aerodynamic  configuration  and  engine  characteristics,  Eq. 
(87)  provides  the  constant  optimum  Mach  number,  say  M «■  0.85.  Then 
0^  (M)  and  K(M)  can  be  evaluated  and  consequently  also  the  value 


o 


i 
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Emax,Say  Emax  “ 8'^*  This  Provides  the  operating  point  on  the  right 
plot  of  Fig.  7.  For  the  example  aircraft  we  can  evaluate  the  in- 
itial wing  loading  0^  by  Eq.  (95).  Assume  that  we  have  the  value 
= 5.  This  translates  into  the  initial  altitude  Z±  * 40,000  ft 
for  cruising  flight.  The  final  altitude,  which  is  higher,  is  ob- 
tained along  the  curve  0^. 

Maximum  Endurance 

In  the  problem  of  maximum  endurance  with  free  range,  we  have 
C2  “ 0.  Then,  the  Eqs.  (76),  (77)  and  (78)  become 


CD  [2  + n(1-n)/2nMKc  + (CD  Kc)m]  = ^-[2  + q(1"n)/2nMK  -(KK  ) ] 
o o M 


(n+l)/2n  M PM  , u>2 

n C1  ~ (CD  + K ~4} 

o M 


Hence,  we  have 


u = M 


rC_ 

2 1 / o 


(96) 


and 


2 

w = M v — - 


A1/2 
K AE 


(97) 


where 


2 + 6Cn"1)/2nMk  + (Cn  k )M 
c D c M 

A * 2 

E 2 + jOi-l)/21^  _ CKR  j 

c CM 


(98) 


50 


In  principle,  the  equilibrium  condition  (74)  and  the  optimum 
relation  (96)  and  (97)  constitute  a system  of  three  equations  for 
evaluating  the  optimum  lift  coefficient,  the  optimum  Mach  number 
and  the  optimum  cruise  altitude  for  the  problem  of  maximum  endur- 
ance. We  will  see  later  that  because  the  effect  of  the  altitude 
on  the  endurance  is  of  the  second  order  a certain  indetermination 
occurs  when  we  evaluate  the  optimum  altitude. 

By  combining  the  Eqs.  (74)  and  (96)  we  have  the  optimum  lift 


coefficient 


This  gives 


(100) 


"D  2/KC, 


Thus,  flight  for  maximum  endurance  is  at  maximum  lift-to-drag 
ratio  to  the  accuracy  provided  by  the  assumption  of  equilibrium 
cruise.  In  terms  of  the  aerodynamic  control,  optimum  cruise  for 
maximum  endurance  is  nearly  the  same  as  in  the  case  of  level  flight. 

Eliminating  to  between  the  Eqs.  (96)  and  (97)  we  have  the  equa- 
tion for  evaluating  the  optimum  Mach  number 


(101) 


Explicitly,  we  have 


(r5-)  -° 

\ max/M 


(102) 
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It  is  seen  from  this  equation  that  optimum  Mach  number  is  con- 
stant. We  have  the  same  remarks  as  given  above  for  the  case  of  max- 


imum range  in  the  stratosphere.  Once,  the  Mach  number  has  been  ob- 
tained from  Eq.  (102),  the  optimum  flight  altitude  is  given  by  Eq. 
(96).  We  can  write  this  equation 


As  before,  with  the  case  of  maximum  range,  we  consider  the 
plots 


(103) 


where 


(104) 


fi  = 


W 


kp . SC 


D 


(105) 


o 

Therefore,  Fig.  7 can  be  used  with  this  new  definition  of  fl 
which  is  the  same  as  the  case  of  maximum  range  in  the  stratosphere. 

This  doesn't  mean  that,  for  a given  initial  weight,  flight  for 
maximum  range  and  flight  for  maximum  endurance  are  performed  at 
the  same  altitude.  In  general,  flight  for  maximum  range  is  at  higher 
Mach  number  and  higher  altitude  than  flight  for  maximum  endurance. 

While  for  the  case  of  maximum  range,  once  the  aerodynamics  and 
engine  characteristics  are  given,  one  can  always  evaluate  the  optimum 
Mach  number  and  the  optimum  altitude,  there  is  a certain  indeter- 
mination in  the  case  of  maximum  endurance  as  has  been  observed  by 
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Miller  and  Koch  (Ref.  1). 

To  explain  this  ambiguity,  let  us  consider  the  optimization 
problem  from  the  point  of  view  of  point  performance.  The  usual 
assumption  is  that  of  unaccelerated  flight.  Hence,  from  the  equation 
in  M in  the  system  (73) 


m = m2(cd  + KC2)  = m2cd 

O 


(106) 


Using  the  equilibrium  condition  (74)  in  this  equation  for  evaluating 
t and  then  substituting  into  the  equation  for  w,  we  have 


5e  <Log  w)  ‘ ' <S7V 


(107) 


At  each  instant,  to  maximize  instantaneously  the  endurance,  we  are 

led  to  minimizing  the  mass  flow,  that  is  the  right  hand  side  of  Eq . 

(107) . This  clearly  shows  the  trade-off  between  the  specific  fuel 

consumption  coefficient  Kc  and  the  lift-to-drag  ratio  C^/Cp  as  has 

been  mentioned  earlier.  Since  this  coupling  is  weak,  we  see  that 

maximum  lift-to-drag  ratio  must  be  used.  Then  the  next  step  is  to 

find  the  minimum  of  K /E  , considered  as  a function  of  the  Mach 

c max 

number,  a condition  expressed  by  Eq.  (102).  Unfortunately,  for  a 

typical  aircraft,  E as  a function  of  the  Mach  number  is  a maximum 
r max 

at  low  Mach  number  where  the  effect  of  compressibility  is  negligible. 
That  is  to  say  Eq.  (102)  is  identically  satisfied  for  a whole  range 
of  low  Mach  numbers.  The  optimum  Mach  number  is  undetermined  and 
subsequently,  from  Eq.  (96)  we  have  a whole  range  of  altitude  for 
optimum  cruise. 
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The  conclusion  is  that  the  effect  of  the  altitude  on  maximum 


endurance  is  not  detected  to  the  first  order  of  approximation.  The 
only  requirement  in  aerodynamic  control  is  that  the  flight  be  per- 
formed at  maximum  lift-to-drag  ratio.  For  each  selected  altitude, 
there  will  be  an  optimum  Mach  number  as  given  by  Eq.  (96). 

To  detect  the  second  order  effect  of  the  altitude,  singular 
variable  thrust  control  is  required.  As  has  been  shown  in  the  pre- 
ceding section  the  varying  thrust  magnitude  is  given  by  Eq.  (57) 
which  contains  the  influence  of  the  altitude  through  the  pressure 
ratio  6.  Since  maximum  endurance  cruise  is  performed  at  low  Mach 
number,  we  can  neglect  the  effect  of  compressibility  and,  from  Eq. 
(98)  take  A = 1,  and  then  have  the  expression  (57)  in  the  form 


T , 2 

kpS  2 + 6(n"1)/2n 


Mk 

c 


•>  CD  ”2 
o 


Using  the  optimum  condition  (96)  we  have 


(108) 


w E 0 . x(n-l)/2n 

max  2+6  Mk 


(109) 


Therefore,  we  can  write  the  equation  in  w in  system  (73) 


h (Log  w)  = " -r~  ( 


„ (n-l)/2n  „ 

max  2+6  Mk 


-) 


(110) 


As  compared  to  Eq.  (107)  we  now  have  the  effect  of  the  altitude 

through  the  pressure  6.  To  minimize  the  mass  flow,  we  are  led  to 

maximizing  the  denominator  on  the  right  hand  side  of  Eq.  (110). 

When  E is  constant,  that  is  when  the  optimum  Mach  number  is  in 
max 


■ I 


the  low  subsonic  range,  we  maximize  the  term  6 ^ Mk^  which, 
because  of  the  optimum  relation  (96) , can  be  expressed  as 


„(n-l)/2n  2w  ,1/2  kc 

6 Mk  = ( -)  t-/,- 

C kP*SC*  61/2n 


(111) 


Maximizing  this  term  is  equivalent  to  taking  the  smallest  value  of 
the  pressure  ratio  6.  Hence  the  optimum  altitude  is  the  highest  one 
possible  provided  that  the  corresponding  Mach  number  as  obtained 
from  Eq.  (96)  is  less  than  the  critical  Mach  number  M^  beyond  which 
the  effect  of  compressible  fluid  is  being  felt.  This  effect  of  the 

* 1/2 

altitude,  at  low  subsonic  range,  is  small  since  the  term  (2W/kp^SCL) 

_3 

is  of  the  order  of  unity  while  a typical  value  of  k^  is  8.7749x10 

_3 

for  a transport  aircraft  (KC-135A)  and  8.9838x10  for  a typical 

fighter  aircraft  (F-4C).  Hence  k /6^2n  varies  from  0.547  k at  sea 

c c 

level  to  k^  at  the  tropopause  and  the  resulting  improvement  in  fuel 
saving  is  small.  On  the  other  hand,  as  the  selected  cruise  altitude 
is  further  increased,  the  optimum  Mach  number  is  greater  than  M^, 
the  effect  of  compressible  fluid  drastically  decreases  and  the 

right  hand  side  of  Eq.  (110)  increases  in  absolute  value  with  the 
altitude.  Hence,  the  mass  flow  increases  and  the  endurance  de- 
creases . 

In  summary,  the  endurance  increases  slightly  with  the  altitude, 
until  a certain  critical  altitude  is  reached.  This  altitude  is  ob- 
tained from  Eq.  (96)  by  using  M = M^  where  M^  is  the  Mach  number  be- 
yond which  the  effect  of  compressible  fluid  is  no  longer  negligible. 


Beyond  this  critical  altitude,  endurance  decreases  as  the  altitude 
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increases. 


Applications 

As  applications  of  the  analysis  we  shall  consider  two  typical 
military  aircrafts,  the  KC-135A  and  the  F-4C.  For  each  aircraft, 
drag  polars  for  a range  of  Mach  numbers  of  interest  have  been  ob- 
tained from  a correlation  of  wind  tunnel  tests  and  flight  test  re- 
sults. For  each  Mach  number,  the  point  of  maximum  lift-to-drag 
ratio  on  the  polar  is  identified.  If  we  require  that  the  theoretical 


parabolic  drag  polar 


CD=CD  +KCL 
o 


matches  the  same  values  for  CT  and  E , then 

L max 

* 

% -r 

o 


(112) 


(113) 


(114) 


For  the  KC-135A  aircraft,  and  K remain  nearly  constant  for  M < 

o 

0.7  and  then  increase  with  the  Mach  number.  The  slopes  of  the  two 
curves  are  positive  and  become  large  as  the  Mach  number  approaches 
unity.  The  variations  are  shown  in  Figs,  8 and  9 with  the  circles 


representing  the  numerical  data  points. 

It  is  always  possible  to  model  the  functions  and  K as  poly- 

o 

nomlals  in  M.  The  inconvenience  of  this  mode  of  representation  is 
that,  in  the  pertinent  formulas  we  have  derived,  the  derivatives 
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dCD  /dM  and  dK/dM  also  need  to  be  represented  accurately.  These  de- 
o 

rivatives  are  positive  and  increase  with  the  Mach  number.  Polynomial 

representation  of  Cp  and  K may  give  negative  slopes  within  certain 

o 

small  intervals.  Let  M^  be  the  Mach  number  beyond  which  we  have  the 

effect  of  compressibility  and  Cp  and  be  the  corresponding  values 

for  Cp  and  K.  Then,  for  the  type  of  curves  such  as  represented  by 
o 

Cp  and  K in  Figs.  8 and  9,  the  following  representations  give  accur- 
o 

ately  the  functions  and  their  first  derivatives  for  M ^ M^ 


cD  = c + q 

o 1 i=l(b.-M) 


(115) 


K = K.  + (M-M  ) Z j-  (116) 

1-1 (e^M) 

The  rationale  of  this  mode  of  representation  is  that  first, 

when  M = M^ , Cp  = Cp  , K = Then  next  when  M approaches  the 

o 1 

smallest  of  the  b^^  and  e^  Cp  and  K increase  rapidly  provided  that 

o 

c^  and  f^  are  positive.  The  number  of  data  points  that  can  be  used 

to  fit  the  curve  Cp  or  K,  besides  the  initial  point  M = M^,  is  3n. 
o 

For  example,  in  the  case  of  the  KC-135A  aircraft,  besides  the  point 
M = 0.70,  we  take  n=l,  and  choose  the  points  M = 0.80,  0.84  and  0.86. 
Then,  each  of  the  two  equations  (115)  and  (116)  provides  three  equa- 
tions for  evaluating  three  constant  unknowns.  The  results  are 


0.0123  + 


2.29x10  (M-0.7) 


(0.8665  - M) 


0.69339 


(117) 


K = 0.05864  + 


7 . 30x10  (M-0.7) 


(0.9742  - M) 


2.884012 


(118) 
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If  the  last  point  is  selected  at  M = 0.90,  the  general  behavior 


of  the  two  curves  for  and  K are  correct  up  to  that  value  of  the 

o 

Mach  number  but  the  accuracy  in  the  low  Mach  number  range  will  de- 


crease unless  we  take  more  data  points  by  using  n=2. 


As  shown  in  Figures  8 and  9,  with  n=l,  for  the  KC-135A,  in  the 


range  of  Mach  number  of  interest  between  0.70  and  0.86  the  Eqs. 


(117)  and  (118)  generate  accurate  numerical  data  for  and  K and 

o 

their  derivatives.  For  M < 0.7  the  values  CL  and  K.  are  used  for 

1 1 

CL  and  K.  Figure  10  plots  the  function  E versus  M with  excellent 
D & r max 

o 

agreement  with  the  data  from  flight  tests.  It  is  seen  that  the  max- 


imum of  E occurs  when  M < 0.70.  Optimum  Mach  number  for  maximum 
max 

endurance  is  selected  in  this  range.  The  optimum  Mach  number  for 


maximum  range  is  given  by  Eq.  (87)  for  flight  in  the  stratosphere. 


Since  k , as  a function  of  the  Mach  number  is  nearly  constant,  the 
c 

condition  is  the  same  as  the  condition  for  the  maximum  of  the  func- 


tion M E which  is  plotted  in  Fig.  11.  From  this  figure  it  is 
max 

seen  that  optimum  Mach  number  for  maximum  range  is  M = 0,77.  Numer- 


ical solution  can  be  obtained  by  solving  the  equation 


M 

/icc 


— > = o 

D 


(119) 


o 

where  and  K are  given  by  the  Eqs.  (117)  and  (118),  Once  the 
o 

Mach  number  has  been  obtained  the  optimum  altitude  for  a given  weight 
is  obtained  from  Eq.  (79),  or  graphically  from  Fig.  7. 

For  the  case  of  maximum  endurance,  as  has  been  shown  above, 
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Fig.  10.  Variation  of  E for  the  KC-135A. 

max 
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the  optimum  altitude  is  the  highest  one  possible  provided  that 
M <^M^.  Since  the  fuel  consumption  for  climbing  to  the  cruise  al- 
titude must  be  taken  into  consideration,  and  since  the  influence  of 
the  altitude  on  performance  is  small  provided  that  the  optimum  Mach 
number  is  selected,  the  procedure  to  follow  is  to  first  select  the 
altitude  and  then  deduce  the  optimum  Mach  number  from  Eq.  (96). 

This  equation  is  plotted  in  Figure  12. 

For  the  F-4C  aircraft,  compressibility  effect  is  felt  for 
M > 0.6  but  stays  small  until  M = 0.7.  Hence,  we  shall  take  M^  = 
0.70.  From  flight  test  data,  the  following  representations  for 

0^  and  K are  accurate 
o 


= 0.0154  + 


2.277x10 

(0.9848-M) 


(M-0.7) 

0.706112 


(120) 


K = 0.155203  + 


0.06573(M-0.7) 


(0.983-M) 


0.442726 


(121) 


The  points  selected  for  evaluating  the  formulas  are  the  points  M = 
0.80,  0.90,  0.95. 

For  the  range  of  Mach  number  of  interest,  between  0.7  and  0.95 

the  functions  are  in  good  agreement  with  flight  test  data  as  shown 

in  Figures  13  and  14.  Figure  15  plots  Emax  as  function  of  the  Mach 

number  while  Figure  16  presents  the  variation  of  M E . From  the 
b r max 

figure,  it  is  seen  that  the  optimum  Mach  number  for  maximum  range 
is  M = 0.85. 

Numerical  solution  can  be  obtained  by  solving  Eq.  (119)  with 
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CQ  and  K as  given  by  Eqs.  (120)  and  (121).  Once  the  optimum  Mach 
o 


number  has  been  obtained  the  optimum  altitude  for  maximum  range  is 
given  by  Eq.  (79)  or  graphically  is  obtained  from  Figure  7. 

Finally,  the  function  (96)  for  determining  the  optimum  Mach 
number  for  maximum  endurance  is  plotted  in  Figure  17  for  the  F-4C 
aircraft.  For  each  selected  altitude,  and  initial  weight,  we  eval- 
uate w and  then  deduce  the  optimum  Mach  number  from  the  graphs. 


Fig.  1 3.  Variation  of  C for  a Typical  Fighter  Aircraft  (F-4C). 
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Variation  of  ME  for  the  F-4C. 
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Fig.  17.  Optimum  Mach  Number  for  Maximum  Endurance  as 
Function  of  the  Altitude  (F-4C). 


SECTION  VII 


COMPLETE  ANALYSIS  OF  THE  CRUISE  PERFORMANCE  PROBLEM 

In  Section  VI,  we  have  used  a rather  engineering  approach  to 
derive  the  optimum  lift  and  thrust  controls  for  the  cruise  perfor- 
mance problem.  The  point  which  can  raise  objections  is  whether  or 
not  it  is  legitimate  to  use  the  altitude  as  a control.  The  arguments 
used  are  based  on  physical  considerations  rather  than  mathematical 
justifications.  On  the  other  hand,  the  advantage  of  the  approach, 
if  it  is  correct,  is  that  we  obtain  explicit  and  simple  law  for  the 
optimum  lift  control  regardless  of  the  end  points  while,  as  usually 
is  the  case  with  an  optimal  control  problem,  the  control  law  can  be 
obtained  only  after  a two-point  boundary  value  problem  has  been 
solved . 

In  this  section,  we  will  corroborate  the  results  obtained  by 
the  simple  analysis  given  in  the  previous  section  by  considering 
the  full  set  of  equations. 

Variational  Formulation 

The  problem  is  formulated  as  an  optimal  control  problem.  The 
state  variables  are  the  dimensionless  quantities  x,  n,  M,  Y and  w. 

The  initial  values  of  the  state  variables  are  prescribed  at  the  in- 
itial time  0=0.  At  the  final  time,  0 = 0^,  the  final  values  of 
the  state  variables  are  partially  prescribed.  It  is  proposed  to 
find  the  optimal  control  law  for  the  lift  coefficient  and  the  di- 
mensionless thrust  t to  maximize  xf  or  0^. 
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"1 


First,  we  have  the  state  equations 


dx  (l-n)/2nM 

de  n * 


V<n+1)/2nMY 


(n-1) / 2n  , , ~ 

^ = 15 [t  - n-1H2(C  + KC2)] 

d0  w D L 

o 


(n-l)/2n  n 9 

a ■ <n  m cl  - “> 


dw 

^ = - K T 
d0  c 


(122) 


- 

dO 


The  Hamiltonian  of  the  system  is 


H = Pc 


+ + KBPn„<"+1>/2"H,  - n-<"+1)/2”(CD  +KC2) 


(,n-l)/2n 

♦ (n’Vc^-w)  + J[n(n-1)/2"PM  - »KCPV) 


o 

(123) 


The  control  is  the  lift  coefficient  and  the  dimensionless  thrust 
t subject  to  the  constraints 


and 


0 1 CL  1 CL  (M> 

max 


0 < t < t (M) 
— — max 


(124) 


(125) 
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The  adjoint  components  are  governed  by  the  adjoint  system 


d_h 

d0 


dp 

r> 

de 


= 0 


= 0 


2nn  ^ - (n-l)n(1-")/2\M  - (n«)K]sPn„<"+1)/2nMY  - 
(n+l)n_(n+1)/2n  — P-  (C  + KC2)  + 

W U L 

O 

-(n+l)/2n  2 

— p [(n+l)M  CL  + (n-l)nw]  - (n-l)n 


(126) 


(n-l)/2n 

w 


„ dpM  (l-n)/2nw  „ (n+l)/2n,.  n(°  1)/2n  , -1„  , w s 

M d r = -pxn  * ' Vnn  ^Y ~ Mpy(ri  CL  + ”2> 

M 

-(n+l)/2n  _ „ 

+ 11 w M PM[CD  (2+CD  )+KCL(2+Km>1  + KcKc  PwT 

o oM  M 

dPY  v (n+1)  / 2n., 

Yd e = -KBpnn  ^ 


P„„-(n+l)/2n 


(n-l)/2n 


^ -(n+1)/2n(C  + kc2}  + Urn  -_V  MC  + n_ 

de  w D L w L w rM 

o 


The  problem  is  solved  by  integrating  the  state  and  adjoint  equa- 
tions, subject  to  the  end  conditions,  while  selecting  at  each  instant 
0 the  controls  and  x such  that  the  Hamiltonian  H is  an  absolute 
maximum. 

With  respect  to  C H is  a maximum  either  for  C = C or 

Lt  L Li 

max 

• . ,<r  when  3H/3CL  - 0.  Explicitly,  we  have 
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With  respect  to  the  linear  thrust  control  t,  we  define  the 


switching  function 


a (n-l)/2n 
6 = n pu  - wK  p 

rM  crw 


Then,  when: 


(128) 


0 > 0,  we  choose  t = r 

max 

0 < 0,  we  choose  t = 0 

0=0,  for  a finite  time  interval,  we  have 
t = variable. 


The  variable  thrust  control  is  determined  by  singular  control  theory. 
We  shall  assume  that  the  fuel  weight  (W^  - W^.)  is  large  enough  so 
that  we  can  use  a suboptimum  climb  to  cruise  altitude  and  consider 
the  flight  as  initiating  with  a sustaining  arc  along  which  r is  var- 
iable. Along  this  arc,  we  have  constantly  © =0,  that  is 


(n-1) /2n  n 
n PM 


wK  p 
c w 


(129) 


Integrals  of  Motion 

The  problem  has  a number  of  integrals 


pe  = ci 


p2  - c2 


(130) 


The  last  equation  in  (122)  is  introduced  so  that  we  can  write  H = 0, 


that  is  with  the  control  law  (127)  , and  0=0 


(n+l)/2n  1/n  (n+l)/n 

c^n  + c2n  m + KBPnn  My 


M2p 


M 


, 2KojC 

w <CD  -KCL+TT-> 
o M 


(131) 


This  relation  is  valid  along  a coasting  arc,  t = 0,  and  along  a var- 
iable thrust  arc,  0 =0. 

At  this  point,  it  is  useful  to  summarize  the  problem.  There 
are  five  state  variables,  n,  M,  y,  w.  At  the  initial  time,  we 
have 


0 = 0,  xi  = 0,  n± , M± , Yi,  wA 
At  the  final  time  we  have 


0^,  x^,  , Mj , Yf , wf 

For  maximum  range  problem,  we  maximize  x ^ while  allowing  0^  to  be 
free.  This  leads  to  Cj  = 0 by  the  transversality  condition.  For  max- 
imum endurance  problem,  we  maximize  0^  while  allowing  x^  to  be  free. 
This  leads  to  s 0 by  the  transversality  condition. 

The  system  of  state  equations  (122)  is  non-linear.  Hence,  op- 
timum trajectory  can  only  be  obtained  by  numerical  integration  using 
some  specified  aerodynamic  and  engine  data.  The  best  we  can  hope 
for  is  to  obtain,  by  some  way,  the  optimum  lift  coefficient  and 
the  optimum  dimensionless  thrust  x as  function  of  the  state  variables 
so  that  upon  substituting  into  the  Eqs.  (122),  optimum  trajectory 
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can  be  derived  by  numerical  integration  of  the  system. 


In  this  respect,  by  Eq.  (127),  the  lift  coefficient  can  be  ex- 
pressed in  terms  of  the  state  variables  if  p and  Mp„  can  be  obtained 

Y M 

as  functions  of  these  variables.  On  the  other  hand,  along  the 
coasting  arc  and  along  the  sustaining  arc,  by  the  Hamiltonian  inte- 
gral (131) , can  be  expressed  in  terms  of  the  state  variables  if 
Mp^  and  np^  can  be  expressed  in  terms  of  these  variables. 

Along  the  sustaining  arc  we  have  the  additional  relation  (129) 
which,  unfortunately  does  not  contain  either  p^  or  p^.  Nevertheless, 
since  the  equation  is  valid  over  a finite  time  interval,  which  is 
the  time  interval  for  sustaining  flight,  we  can  take  its  derivative 
to  have  an  additional  relation.  It  is  known  that  the  linear  control 
does  not  appear  with  this  first  derivative  (Reference  4).  By  taking 
the  derivative  of  Eq.  (129)  with  respect  to  0,  and  neglecting  the 
effect  of  T),  assumption  which  is  rigorously  exact  for  n=l,  we  have 

M^p 

c2n1/n  M + KBPnn(n+1)/nMY  = {CD  [2  + n(1_n)/2nMKc  + (CD  Kc)m1 


KC2[-^-+  n(1-n)/2nMK  - (KK  ) ]} 

L m2^  c c M 


(132) 


The  three  relations  (129),  (131)  and  (132)  among  wPwi  PP^  and  Mp^ 
allow  the  elimination  of  two  of  these  variables.  We  need  to  find 
two  more  integrals,  at  least  one  of  them  involving  yp^.  In  this 
respect,  the  non  integrable  equations  of  the  system  (126),  with 


2rv 


d(npn> 


d0 


. (l-n)/2n  , . (n+l)/2n 

(,n-l)n  C M + (n-l)K  p n MY  - 

i Bn 


rnm,r(n+1)/2n  M ?M  rr  irr2  2(n-l)  wKCL. 
(n+l)n  — [cD  - kcl  - -^y  ~2~ 1 

o M 


(n-l)n(n_1)/2n  ^ T 
w 


- -C2nW-"),2nM  - KBPnn<,,'-1>/2"M,  + 


-tn+l)/2n  2KwC 

„ M PM[CD  (1+CD  )-KCL(1"KM)  “ ~1~]  + (133) 

° °M  M 


(,n-l)/2n 


Pm(1+Kc  )t 


w c, 


M 


d (yp  ) , ,lWo  o -(n+l)/2n  0 0 wKCT 

V • ^ ^ 


d(wpw)  -(n+l)/2n  2 2 

~df~ " ” “ — w — mVcd  ~kcl> 
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These  equations  are  valid  for  sustaining  flight,  and  also  for 
coasting  flight  with  t = 0. 

The  only  case  where  we  can  obtain  an  exact  integral  of  this 
system  is  the  case  for  flight  in  the  stratosphere.  With  n=l,  we 
have  the  equations  in  np^  and  wp^ 


d(npn) 

x2p„ 

2 

d 0 

w 

CD 

0 

d(wpw) 

h2P„ 

<s  -kc^> 

r* 

d0 

V 

(134) 
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Therefore 


l 


d(wp  ) d(np  ) 

- = 

de  de 

and  upon  integrating 


wp  = np  + C0 
*w  r 0 3 


(135) 


(136) 


where  is  a new  constant  of  integration 

When  n^l,  in  order  to  obtain  the  required  additional  integrals 
some  simplifying  assumptions  have  to  be  made.  If  the  optimum  re- 
lation (132)  for  sustaining  flight  is  used  to  eliminate  the  and 
terms  in  the  equation  in  Mp^  in  system  (133),  we  have 


d(MpM) 

de 


mVm 

0) 


<CD  - KCJ)  + 2 
o 


-(n+1) / 2n 


)tnT-M2(CD  +KC2)] 
o 


(137) 


Now  if  we  assume  that  the  Mach  number  stays  nearly  constant, 
that  is  at  each  instant  of  flight  dM/d6  ~ 0,  we  have  the  assumption 
of  quasi-unaccelerated  flight 


DT  = M2(Cq  + KC2)  (138) 

o 

Using  this  assumption,  not  as  additional  relation  for  providing 
the  thrust  law,  but  only  as  a simplifying  device  to  obtain  approx- 
imate integrals  of  the  motion,  the  Eq.  (137)  is  reduced  to 


d(MpM) 

de 


0) 


(139) 
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Replacing  Kc  by  its  expression  obtained  from  the  singular  re- 


lation (129) 


d(MpM) 


(n-l)/'2n  M3Pm2  , 

n — (c  - kc2) 

0)  wp  D L 

w o 


(140) 


Combining  this  equation  with  the  equation  for  wp^  in  system 
(133),  we  have 


d(wpw)  d(MpM) 


(141) 


The  integration  yields 


wp  = Mp  + C. 
w M 4 


(142) 


where  C^  is  another  constant  of  integration. 


In  system  (133) , by  eliminating  the  K_  terms  between  the  two 

D 


equations  for  np^  and  yp^  we  have 


d(np  ) 


+ (n-1) 


d(yp  ) 


(n-l)C2n(1_n)/2nM  - 


J -(n+l)/2n  1 i \ n ” ( n"f“l ) / 2 n 

^ M2p  (C_  - KC2)  - 

w rM  D L w 

o 

Pm[ht  - M2(Cd  + KC2)] 
o 

Combining  this  equation  with  the  equation  for  wp^,  we  have  the  exact 


differential  relation 
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d("<>»)  , (.-1)  ln-1)  . (l-n)/2n.. 

d0  ’ d0  2n  d0  2n  ^2n 


“(n+l)/2n  7 7 

2n* Pn1"'  - M (CD  + KCL)] 

o 


(146) 


If  n = 1,  we  can  integrate  the  equation  and  have  the  exact  integral 
(136).  When  n ^ 1,  we  can  use  the  condition  of  unaccelerated  flight 
(138)  to  reduce  the  equation  to 


d(nPn)  d(“1’»)  , (n-1)  d<YPv)  (n-1)  . (l-n)/2n,, 

d0  d0  2n  d0  2n  ^2n 


(144) 


The  equation  can  be  integrated  in  the  case  of  maximum  endurance 
with  free  range,  C 2 = 0.  This  gives 

(n-1) 


wp„  = np„  + 


2n 


+ C, 


(145) 


We  need  one  more  integral  for  the  case  of  maximum  range  with 
free  final  time,  = 0.  This  integral  is  not  necessary  to  obtain 
the  lift  control  law  if  we  only  consider  the  sustaining  trajectory 
since  when  = 0 we  can  combine  the  two  relations  (131)  and  (132) 
to  have 


Ml  - n(1-n>/2Ve  + (KKc)„]C2  - ^ CL  + 

M 

CD  [1  + n(1~n)/2nMKc  + (CD  kc)m)  = 0 


(146) 


This  quadratic  equation  can  be  solved  for  the  optimum  lift  co- 
efficient CL  in  terms  of  the  state  variables  n,  M and  w. 
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The  constants  of  integration  must  be  evaluated  by  using  the 
boundary  conditions.  This  requires  a complete  integration  of  the 
state  and  adjoint  equations  over  the  entire  trajectory  which  may  in- 
clude a maximum  thrust  arc.  An  analytical  integration  of  these 
equations  is  obviously  not  possible  in  the  present  case.  Neverthe- 
less, the  integrals  obtained  are  valuable  because  of  the  following 
reasons : 

1.  They  reduce  the  order  of  the  system.  Each  integral  allows 
the  elimination  of  one  differential  equation. 

2.  They  may  provide  relations  among  the  state  variables  or 
relations  among  the  state  and  control  variables.  In  the  former  case 
these  relationships  show  explicityly  that  the  optimum  trajectory 
must  lie  on  a certain  surface.  For  example,  in  Section  V,  it  is 
shown  that  optimum  singular  arc  is  the  curve  as  given  by  Eq.  (54). 

In  the  later  case,  the  relations  obtained  give  the  optimum  control 

as  function  of  the  state  variables.  This  is  illustrated  by  Eq.  (146) 
above  for  the  lift  coefficient  in  the  case  of  maximum  range  and 
later  on  in  this  section  in  Eq.  (212)  for  the  lift  coefficient  in 
the  case  of  maximum  endurance  and  in  Eqs.  (199)-(201)  and  (221)-(223) 
for  the  thrust  control, respectively, in  these  cases. 

3.  When  used  with  an  independent  numerical  program,  the 
exact  integrals,  if  they  exist,  can  be  used  to  verify  the  accuracy 
of  the  numerical  scheme. 

4.  Finally,  in  some  favorable  cases,  the  integrals  obtained, 
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although  not  sufficient  in  number,  may  provide  approximate  optimum 
control  law  which  can  be  used  as  a first  estimate  for  an  iteration 
procedure. 

As  illustration  of  this  last  statement  we  shall  now  use  the  in- 
tegrals obtained  to  corroborate  the  results  of  the  preceding  section. 
First,  let  us  consider  the  case  of  flight  in  the  stratosphere. 

Let 


A = Mp, 


M 


B = wp 


w 


(147) 


c = np. 


D = YP. 


where  here  A is  different  from  the  A defined  above  for  AR  and  A^. 

For  flight  in  the  stratosphere,  n=l,  and  we  can  write  the  Eqs.  (129), 
(136)  and  (142)  as 


A = MK  B 
c 


B = C + C, 


B = A + C, 


(148) 


This  linear  system  can  be  solved  for  A,  B,  C in  terms  of  MKc> 
that  is, in  terms  of  the  Mach  number, provided  that  MKc  ^ 1. 

Now,  in  consistence  with  the  assumptions  used  in  Section  VI  we 
assume  that  the  flight  path  angle  y not  only  is  small  but  varies  very 
slowly  along  the  flight  path.  Hence,  we  have  the  equilibrium  con- 
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dition 


0) 


M2 


(149) 


Using  this  condition,  we  rewrite  the  integrals  (131)  and  (132) 
for  n=l 


r+c2  + +¥> 

o M 


(150) 


and 


C2  + W = « {CD  t2  + MKc  + (CD  Kc)m]  - 


(151) 


Kw 


T 12  + MK c " (KKc)M11 

M 


The  term  in  is  small  because  of  the  coefficient  y.  Furthermore, 
it  will  be  shown  later  that  C is  small,  hence  KgCy  can  be  neglected. 

We  can  also  keep  this  term  and  assume  that  y is  constant,  an  assumption 
which  is  consistent  with  the  equilibrium  cruise  condition.  Then, 
since  in  the  Eqs.  (150)  and  (151),  A and  C are  functions  of  the  Mach 
number,  these  equations  consist  of  a system  of  two  equations  for  the 
variables  0)  and  M.  Hence,  for  flight  in  the  stratosphere  both  the 
Mach  number  M and  the  wing  loading  w are  constant.  These  results 
agree  with  those  in  the  previous  section  and  are  valid  for  both  the 
maximum  range  and  the  maximum  endurance  problems.  This  should  be 
noted  here  that  they  are  only  approximate.  It  will  be  shown  later 
in  a more  rigorous  analysis  that  the  Mach  number  slowly  decreases 
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along  the  flight  path.  From  system  (148),  under  the  present  assump- 
tions, it  appears  that  A,  B and  C are  also  constant.  From  any  one 
of  the  equations  for  these  variables,  say  Eqs.  (134),  we  see  that 


C 


L 


(152) 


The  flight  for  both  the  maximum  range  and  the  maximum  endurance  in 
the  stratosphere  is  at  maximum  lift-to-drag  ratio.  From  the  equation 
for  p^  in  system  (126),  we  have,  for  n=l 


- - KbCK  (153) 


On  the  other  hand,  from  the  optimum  relation  (127)  for  the  lift  co- 
efficient 


pT  ■ 2kacl  ■ r4-  <154> 

' max 

The  right  hand  side  of  this  equation  is  function  of  the  Mach  num- 
ber and  therefore  is  approximately  constant.  Hence  p^  is  near  con- 
stant and  from  Eq.  (153)  it  appears  that  the  constant  C = np^  is 
small  and,  as  has  been  mentioned  above,  that  ICgCy  can  be  neglected. 

For  the  case  of  maximum  range  with  free  final  time,  = 0,  com- 
bining the  two  equations  (150)  and  (151)  we  have 


a)  = M 


2 


(155) 


For  the  case  of  maximum  endurance  with  free  range,  C2=0,  using 
the  Eq.  (151)  with  the  term  l^Cy  neglected,  we  have 


0)  = 


(156) 
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These  results  are  in  perfect  agreement  with  those  in  Section 


1 


VI. 


Next,  we  consider  the  case  of  maximum  endurance  in  the  tropo- 
sphere,  n ^ 1.  With  A,  B,  C and  D as  defined  in  (147),  we  write 
the  integrals  (129),  (142)  and  (145)  which  are  valid  in  this  case 


A = n(1"n)/2nMK  B 


B = A + C, 


B = C + ("~1)  D + C 
2n  3 


(157) 


Also,  we  have  the  Eqs.  (131)  and  (132),  written  with  = 0 and  using 
the  equilibrium  cruise  condition  (149) 


C.  v /0  (n-l)/2n  2 

— + K n (l~n)/2n_  

M V y u M D m4  ; 

o M 


(158) 


and 


(1-n)/2nCl  , I*"  “/2"  A{  [2  + „(1-n)/2nMKc  + (CD  K )„] 

o o 


- [2  + n(1-n)/2nMK  - (KK  ) ]} 

„4  c c M 

M 


(159) 


The  optimum  lift  control  law  (127)  is  written  in  terms  of  A and 


D as 


Al_ 


max 


(160) 


If  y is  considered  as  constant,  an  assumption  which  is  in  line  with 
the  equilibrium  cruise  condition,  then  the  Eqs.  (157)-(160)  constitute 
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a system  of  6 equations  for  the  unknowns  A,  B,  C,  D,  M,  a)  .and 
ti (1  Hence,  the  first  six  variables  are  functions,  of  n,  that 

is  of  the  altitude.  This  influence  of  the  altitude  is  small.  For 
example,  let  us  consider  an  altitude  gain  of  about  2,000  ft, which 
is  normally  over  a few  hundred  miles  in  range  and  in  time  interval 
long  enough  to  insure  the  assumptions  of  equilibrium  cruise  and 
unaccelerated  flight.  Then,  we  have  the  ratio 


(l-n)/2n 


/ Pi  \ (n-l)/2n 

Vf ) 


(161) 


With  n = 1.235  and  an  altitude  climb  of  2,000  ft.,  this  ratio  varies 
from  1.008005  at  20,000  ft  to  1.009154  at  36,000  ft.  Hence,  as  a 
first  order  approximation  we  can  neglect  the  effect  of  the  altitude 
and  deduce  that  the  quantities  A,  B,  C,and  D are  quasi-constant. 

If  A or  B are  considered  as  constant,  then  from  their  equations 
we  deduce 


C 


L 


(162) 


A result  which  is  in  agreement  with  the  control  law  derived  in  Sec- 
tion VI.  If  D is  considered  as  constant,  from  the  equation  in 
yp^  = D in  system  (133)  with  the  equilibrium  cruise  condition 
applied,  we  are  led  to  neglecting  the  term  n^^nCMY.  Now,  if 

C = np  is  considered  as  constant,  then  from  the  first  of  the  Eqs. 
(133),  using  the  condition  for  unaccelerated  flight  (138)  and 
equilibruim  cruise  (149) , we  have 
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, (l-n)/2n  2nn(n"1)/2n  „„2 

(n-DKgn  Cy A(Cd  - KCL) 


(163) 


From  the  fact  that  the  left  hand  side  is  small,  we  have  the  control 
law  (162). 

Finally,  we  consider  the  case  of  maximum  range  with  free  final 
time,  = 0,  for  flight  in  the  troposphere,  n ? 1.  As  mentioned 
above  there  is  one  more  integral  to  be  found  for  this  case.  If  the 
approximate  condition  (149)  is  used  in  the  exact  relation  (146) , we 


KC2  1 + n<1~n)/ZnMKc  + (Cp  Kc)„ 
L _ 2 

CD  3 + n(1_n)/2nMK  - (KK  ) 

o c c M 


(164) 


Using  the  integral  (131)  with  = 0 to  eliminate  the  terms  C 2 and 
Kb  in  the  equation  for  np^  in  system  (133)  and  then  applying  the 
assumptions  (138)  and  (149) , we  have 


d(2nnp  ) -(n+l)/2n  , , 

— 10^“  = " " w M PM[(n+1)CD  " (3-1>KCL]  (165) 

o 


Combining  this  equation  with  the  equation  in  wpw  in  system 


(133) 


f2nnPn  - (n+i)wpw]  = ys=pL  n"(n+1)/2nM2PMKC2  (166) 

On  the  other  hand,  by  subtracting  the  two  equations  in  np^  and 

wp  in  system  (133)  using  the  same  assumptions,  we  obtain 
w 
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Again,  if  Eq.  (131)  with  = 0 is  used 

d -(n+l)/2n  ~ 2 

it  [2n<"pn  -”!>„>'■  '"-D  2 7 » P„<CD  + kcl> 

o 

Dividing  the  two  equations  (166)  and  (168)  and  noticing  that 


(167) 


(168) 


KVCd  =AR 

o 


(169) 


we  have 


d[2nnpri  - (n+l)wpw]  2AR 


d[2n(riPn  - WPW)]  1 + \ 


(170) 


As  shown  above,  the  effect  of  the  altitude  through  the  term 
n<l->/2"MK  on  the  right  hand  side  of  this  equation  is  negligible, 
not  only  that  n)/^n  near^y  constant  but  also  because  MKc  is 
negligibly  small.  Hence,  the  right  hand  side  of  the  Eq.  (170)  is  a 
function  of  the  Mach  number  and  if  the  condition  of  unaccelerated 
flight,  M = constant,  is  applied  solely  for  the  purpose  of  perfor- 
ming the  integration,  we  have 


f(3n-l)A  - (n+1) ] 2n(AR  - 1) 

~ 1 -L  a \ np„  + Cc 


(1  + V 


"w  (1  + a^  "‘n  5 


(171) 


where  is  another  constant  of  integration.  This  integral  replaces 
the  Eq.  (145)  for  the  case  of  maximum  range  in  the  troposphere. 

We  now  can  summarize  the  integrals  of  interest  in  this  case, 
that  is , the  integrals  (129),  (142)  ,and  (171)  as 


A = n(1“n)/2nMK  B 


B = A + C, 


(172) 
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[(3n-l)AR  - (n+1)]  2n(AR-l) 

(1+A^)  B “ (1+A^  C + C5 


where  we  recall  that  AR  is  defined  by  Eq.  (164).  This  relation  re- 
sults from  combining  the  two  integrals  (131)  and  (132)  so  that  we 
can  write  it  with  one  of  these  integrals,  say  Eq.  (131)  with  0^=0 


(173) 


(l-n)/2n  (l-n)/2n  _ n (n-i) /2n  K 2 

c2n  + KBn  Cy  - — A (CD  + -*-) 

O M 


(174) 


If  the  flight  path  angle  y is  considered  as  constant,  the  equations 
(172),  (173)  and  (174)  constitute  a system  of  5 equations  for  the 
6 variables  A,  B,  C,  0),  M,and  n)/2n.  Hence,  again  we  have  the 
dependence  of  the  first  five  variables  on  the  altitude  through 
r/l  n^^n.  Since  it  has  been  shown  that  this  term  stays  nearly  con- 
stant, the  variables  A,  B, and  C are  quasi-constant.  If  C = np^  is 
considered  as  constant,  then  from  Eq.  (165),  we  have  the  approximate 
optimum  control  law  for  the  lift  coefficient 


/n+1  v/_^c 

3n-l  K 


/n+]__  * 

3n-l  \ 


(175) 


in  perfect  agreement  with  the  law  derived  in  Section  VI  for  the  case 

of  maximum  range  in  the  troposphere.  If  we  consider  A or  B as  con- 

* 

stant,  then  from  their  equations  we  have  » C^.  The  effect  of  n 
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on  the  optimum  lift  coefficient  as  given  by  Eq,  (175)  is  lost. 


I 


The  rationale  behind  the  selection  of  C as  constant  over  A and 
B is  as  follows.  From  the  first  of  the  Eqs.  (172),  since 
n(1'n>/2X  <<  1,  we  have  A <<  B.  By  the  same  arguments  as  used 
above  for  flight  in  the  stratosphere  we  can  show  that  C is  negligibly 
small.  Then  the  term  n^2nCy  in  Eq.  (174)  is  smaller  and  can 

be  neglected.  From  this  equation,  since  we  maximize  the  range,  we 
can  normalize  the  adjoint  variables  by  taking  = 1.  Then  we  can 
write  the  equation  as 


A 


(l-n)/n 
con 


(176) 


By  this  expression,  A has  a finite  value  not  near  zero.  Therefore 
in  decreasing  order  of  magnitude,  we  have  B >>  A >>  C and  it  is  most 
accurate  to  use  C as  a constant  to  obtain  the  optimum  lift  control 
law. 


In  conclusion,  we  notice  that  the  maximum  lift-to-drag  ratio 
law  and  the  constancy  of  the  Mach  number  M and  the  wing  loading  u> 
are  only  approximate.  The  purpose  of  the  analysis  is  simply  to 
justify  the  engineering  approach  of  Section  VI.  We  shall  next  de- 
rive the  main  results  of  this  report,  namely  obtaining  the  exact 
forms  for  the  lift  control  and  thrust  magnitude  for  cruising  flight 
to  obtain  either  maximum  range  or  maximum  endurance. 

Singular  Aerodynamics  and  Thrust  Magnitude  Controls  for  Maximum  Range 
The  thrust  magnitude,  as  given  by  the  condition  (138)  of  un- 
accelerated flight,  is  only  an  approximation  for  the  purpose  of 
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simplifying  the  adjoint  equations  to  make  them  integrable.  To  have 
the  exact  expression  for  the  variable  optimum  thrust  control,  we 
start  with  the  singular  condition  (129)  rewritten  here  for  conven- 
ience 


(n-l)/2n 
n p 


M 


K wp 
c rw 


(177) 


This  condition  is  identically  satisfied  along  the  main  portion 
of  the  trajectory,  namely  the  singular  sustaining  arc.  By  taking 
its  derivative  with  respect  to  0 and  neglecting  the  effect  of  o we 
have  obtained  the  condition  (132).  If  we  retain  the  effect  of  n» 
then  by  taking  the  derivative  of  Eq.  (177)  using  the  state  equations 
(122),  the  adjoint  equations  (126),  the  optimum  lift  condition  (127) 
and  the  singular  condition  (177)  itself,  we  have  the  following  re- 
lation 


c2  + VV  + -T?  KBMpMY  ■ " 


-1/n 


-Hp^Co  [2+nH-’)/2Vc+(CD  Kc>„] 

o o 


- ^[4^  + n(1_n)/2nMK  - (KK  ) ]} 

L i C CM 

M CL 


(178) 


It  is  known  that  the  linear  thrust  control  does  not  appear  in  this 
derivative.  Now,  by  replacing  C by  its  optimum  expression  (127)  in 

Li 

this  equation  and  then  by  taking  the  derivative  we  will  have  a new  re- 
lation in  which  the  thrust  magnitude  x will  appear  linearly.  It 
should  be  noted  that  the  optimum  singular  thrust  control  obtained 
this  way  is  exact  in  the  sense  that  the  process  does  not  involve 
any  approximation  except  the  approximation  of  small  flight  path  angle 


1 


and  results  in  the  state  equations  (122),  which  we  consider  as  exact 
equations  forming  the  basis  for  the  present  analysis.  It  is  only 
after  the  expression  for  the  thrust  has  been  obtained  that  we  use 
the  equilibrium  cruise  condition  to  simplify  it  for  practical  use. 
The  process  of  taking  the  derivative  of  Eq.  (178)  is  tedious  but 
straightforward.  Here  we  shall  only  consider  the  case  of  flight  in 
the  stratosphere.  With  n = 1,  we  have  the  simplified  expression 
MPv 


S + W ' -ir  <S  12  - mc  + <S  Vm>  - fr  cl 

o o M 


KCL[Mkc  - (Kkc)M] } 


(179) 


We  notice  that  K 


k for  n = 1. 
c 


On  the  other  hand,  we  have  the  Hamiltonian  integral  (131)  for 


n = 1 


^1  MPM  y 

— + C„  + K T)P  Y = — (C_  - KC.  + -z-  C ) 

M 2 B n cu  D LW2L 

o M 


(180) 


For  the  case  of  maximum  range  with  free  final  time,  = 0,  by  com- 
bining the  two  equations  (179)  and  (180)  we  have  the  quadratic  equa- 
tion for  C. 


K[1  + «k  > - »c)  <£  - X Cl  + S II  + «i>  T>c>„  + «kc)  ' 0 

M o o 


(181) 


This  equation,  upon  solving,  gives  the  explicit  expression  for 
the  optimum  lift  coefficient  for  maximum  range  in  the  stratosphere. 
Let 
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K = 1 + (Kk  )u  - Mk 
1 c M c 


K2  ' 1 + <CD  kc>M  + ““c 
o 


(182) 


and  write  the  equation  for  the  optimum  lift  coefficient 


- rf  CL  + CD  K2  ■ 0 
M o 


(183) 


Solving  for  C 


2u) 


K1K2CD 


"L  2 1 — 

L Kit 


[1  ± VI  - 


-2.  «4 


M4] 


(184) 


4Ko) 


There  is  some  inconvenience  in  using  this  formula.  First, 
there  is  the  ambiguity  in  the  (+)  sign.  If  the  singular  arc  starts 
at  a point  along  a climb  to  the  optimum  cruise  altitude,  then  the 
initial  value  for  is  large  and  the  (+)  sign  must  be  used.  But 
each  time  the  quantity  under  the  square  root  passes  through  zero, 
we  have  to  change  the  sign  in  front  of  the  radical.  Next,  as  can 
be  seen  in  the  optimum  control  law  (127),  the  lift  coefficient  CL 
is  proportional  to  the  adjoint  component  p^.  Experience  in  numerical 
optimization  has  shown  that  this  variable  is  nearly  constant  along 
the  main  portion  of  the  flight  trajectory,  namely  along  the  singular 
arc,  . :'.cr  he  thrust  has  a nearly  constant  intermediate  level,  but 

undergoes  drastic  variations  at  the  two  end-points.  This  can  be  in- 
terpreted as  that  the  lift  coefficient  C varies  nonlinearly  at  the 

L 

two  end-points  in  order  to  satisfy  the  Initial  and  the  final  condi- 
tions but  stay  nearly  constant  along  the  main  portion  of  the  trajec- 
tory. For  this  type  of  behavior  of  the  variable  p_^,  the  initial 
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value  p , that  is  the  initial  value  C , is  extremely  sensitive  to 
Yi  1 

the  subsequent  trajectory  generated.  In  other  words,  if  the  formula 

(184)  is  to  be  used,  the  exact  initial  point  u>^,  where  the  sing- 
ular arc  is  initiated  must  be  accurately  evaluated  in  order  to  pro- 
vide accurate  initial  value  for  C^.  Any  small  deviation  from  the 

exact  initial  lift  coefficient  C can  result  in  steep  climb  or 

Li 

steep  descent  and  hence  not  lead  to  the  required  final  condition, 

As  will  be  seen  later,  the  usefulness  of  the  explicit  equation 
(183)  for  the  optimum  lift  coefficient  CL  is  that  it  is  used  to  gen- 
erate the  optimum  singular  thrust  control  t.  As  for  the  optimum 
lift  coefficient  itself,  for  all  practical  purposes,  it  is  best  to 


use  the  approximate  control  law 


(185) 


Using  this  control  law,  and  the  approximation  of  dy/d9  * 0,  that  is 


“ = clm 


(186) 


in  the  optimum  relation  (183),  we  obtain  the  approximate  equation 
for  evaluating  the  average  optimum  Mach  number  for  cruising  flight 


K1  + K2  = 4 


Explicitly,  from  Eqs.  (182),  we  have 


(kV?m  ^ 2 


(187) 


(188) 


For  a given  aircraft  aerodynamic  and  engine  characteristic,  this 
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equation,  upon  solving,  provides  the  average  optimum  Mach  number  in 
cruise.  The  optimum  altitude  for  cruise  is  obtained  as  function  of 
the  optimum  Mach  number  by  combining  the  two  equations  (185)  and 
(186) 


a)  » M 


(189) 


In  summary,  for  maximum  range,  a first  order  optimum  cruise  sol- 
ution can  be  obtained  as  follows.  First,  Eq.  (188)  is  solved  for  the 
optimum  cruise  Mach  number.  Next  the  optimum  cruise  altitude  is  ob- 
tained from  Eq.  (189).  The  cruise  trajectory  at  this  altitude  using 
constant  optimum  Mach  number  provides  essentially  the  maximum  range. 

Constant  Mach  number  cruise  is  only  an  approximation.  Rigor- 
ously speaking,  the  optimum  Mach  number  and  hence,  the  optimum  thrust 
coefficient,  slowly  decrease  along  the  sustained  flight  path.  To  de- 
rive the  optimum  dimensionless  thrust  t,  we  use  the  optimum  relation 
(127)  to  rewrite  the  equation  (183)  as 


MKKlP^  - 8K2nwpypM  + 4K2K2Cd  M3pM2  = 0 

o 


(190) 


Upon  taking  the  derivative  of  this  equation,  the  thrust  magnitude 
will  appear  linearly.  For  convenience,  we  reproduce  here  the  deriv- 
atives of  the  state  and  adjoint  variables  involved  for  the  case  n * 1 


ie  - 


f . i[„,  - m2(cd  +k<£] 

O 


r 


§ - -kcT  <191> 


de  " - W 

k P 

dP  MP  c rM 

3 r ■ -WV  + iticd  (2  + SJ+  **A  - ^ cl>+  “Sr  ■” 

o °M  M 

Then,  by  taking  the  derivative  of  Eq.(190),  using  the  Eqs.  (191)  and 
the  Hamiltonian  integral  (180)  with  = 0,  with  the  notation  for 
logarithmic  derivative,  we  have 

161t\rf  F 

'“iV1  + + V H + 4K  CD  WO  + 2KM  + CD  + V1 

M O ow  M 

M 

x^tnT  - M2(Cd  + KC2)]  - 8K2 (KgUMy  - - 2KbKM(MK1py  - 4Kcop.rt)nPn  + 

Mp 

8K2<S  ^Pm  - »PJ  *(-^[CD  U + CD  ) + KC2(1  + V ' 2 S.1  + 


V" 


nt)  = 0 


(iy2) 


Since,  by  relation  (127) 


PY  = 2KMpMCL 


we  can  make  tne  substitution,  and  upon  simplifying  by  4K  p„,  we 

M 


(193) 


[KKjM  ll  + ^ + K )C2  - 810^^  + K2CD  M2(3  + 21^  + Cp  + K )] 

M o o M 

M 

1 2 2 2 2 ^Pn 

X^[nx  - M (CD  + KCl)]  - - kcnT)  + 1^(20)  - K^M  + 

o PM 

^(CD  K2M2  - 2toKCL)  x[M2Cd  (1  + CD  ) + KC^M2(1  + 1^)  - 4KwCL  + 


k nx]  = 0 

Si 


(194) 


This  expression  is  exact  and  can  be  solved  for  nx  which  appears 
linearly  in  the  equation. 

For  the  ratio  np  /Mp..,  we  can  use  the  exact  singular  relationship 
q M 


Mp. , = Mk  wp 
M c w 


(195) 


and  the  exact  integral  (136)  for  n=l 


wpw  = npn  + C3 


(196) 


We  also  have  the  Hamiltonian  integral  (180)  with  = 0 


C2  + W = st(cd  - «£  + CL> 


(197) 


Solving  the  system  of  Eqs.  (195)-(197)  for  qp^  and  Mp^  and  forming 
the  ratio 


MpM  * u)Mkc(C2/C3  - I^y)  lC^  " Mkc(CDQ  " KCL  + CL)]  (198) 


Upon  substituting  this  expression  into  Eq.  (194),  and  then  replacing 
CL  by  its  expression  as  given  in  Eq.  (184),  we  can  solve  for  the 
optimum  thrust  control  nr  in  terms  of  the  state  variables  M,  w,  y 
and  one  constant  of  integration  C^/Cy  In  practice,  it  is  accurate 
to  neglect  the  term  np^,  and  also  the  term  K^y  in  Eq.  (194).  Then, 
the  resulting  expression  for  qr  depends  solely  on  the  two  variables 
M and  w, as  is  the  expression  for  the  lift  control  C^.  Explicitly, 
we  have 


where 


(199) 


P = 4u>M2KCl(2  - K + K )(CD  + KC2)  + M4(K2 

' rl  o o M 


+ K2Cd  M4(Kx  + K2)[(Cd  + KcJ)  - 2(Mkc/K2)KC2] 
o o 


-K  )(c  +KcJ) 

M o 

(200) 


and 

Q = KK1M2C2(2  - K + K ) + K2Cd  M2(2  + K2  + K2  ) (201) 

n Uo  M 

2 

In  the  expressions  for  P and  Q,  KC^  can  be  expressed  in  terms 

of  by  using  Eq.  (183),  but  it  is  more  convenient  to  leave  them 

in  the  present  form.  Then,  since  Mkc/K2  is  small,  it  is  seen  that 

2 

r|T  is  proportional  to  CL  = Cn  + KCL  . If  in  P and  Q,  C.  is  replaced 

DDL  Li 

O 

by  its  approximate  optimum  expression  (185)  and  the  small  term  Mkc/  K2 
is  neglected,  we  have  the  following  approximate  expression  for 
the  optimum  singular  thrust  control 
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11T 


8u>  /KC  (2  - Kn  + K,  ) + 2K,C  M (K  + K,  + K,  - K.  ) 

o Hi  o 1 Hi Hi 

Kx(2-K  +K  ) + K„(2  + K.  + K„  ) 

1 >1  M 


(202) 


Furthermore,  if  the  equilibrium  condition  (186)  is  used  as  an  approx- 
imation for  u>  we  have  the  simplified  expression 


4(2  - K.,  + K,  ) + 1^(1^  + K,  + K0  - K,  ) 


nr  = 


M 


Hi  Hi 


K.  (2  - K + K.  ) + K,(2  + Kj  + K„  ) 

M Hi 


(2Cn  M ) (203) 


Finally,  it  has  been  said  that  the  Mach  number  is  not  constant 
along  the  flight  path.  That  is  because  the  thrust  magnitude  n't  as 
given  by  the  Eqs.  (199)-(201)  is  slightly  less  than  the  drag.  Using 
these  equations,  and  the  Eq.  (183)  for  the  lift  coefficient  in  the 
equation  for  M,  with  n = 1,  in  system  (122),  we  have  after  some  man- 
ipulation 


dM 

d6 


2kc%  <K1  + *2>m5kCL 
o 

gjQ 


(204) 


This  equation  is  the  exact  equation  for  M.  It  is  seen  that 

the  Mach  number  decreases  at  a slow  rate  because  k is  small.  To 

c 

have  an  idea  about  the  reduction  of  the  Mach  number,  we  combine  this 
equation  with  the  equation  for  the  weight  to  have 


dw 

dM 


wP 

2Cd  (Kx  + K2)M5Kc£ 


(205) 


o 

If  in  the  expression  (200)  for  P we  neglect  the  term  Mk^  as  com- 
pared to  the  drag  coefficient  and  then  use  the  approximate  solutions 
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T 


T 


(185),  (186)  and  (187),  we  have  the  approximate  equation 


dW  _ bdM 
W " M 


(206) 


where 


+ k2) 


+ K K + K,K„  ] 
1 *M  M 


(207) 


Holding  b constant  for  the  integration  of  Eq.  (206)  we  have  the 
approximate  relationship 


\ /M1/b 

Mf  'l  V 


(208) 


Figure  18  presents  the  plots  of  K^,  K2  and  + K2  for  the 

KC-135A  using  the  functions  for  CQ  and  K as  developed  in  Section  VI. 

o 

The  point  where  + K2  = 4 gives  approximately  the  optimum  Mach  num- 
ber for  maximum  range.  At  that  point  we  have  M = 0.774  and  b = 
8.938594.  Likewise,  Figure  19  presents  the  plots  of  the  same  func- 
tions for  the  F-4C  using  the  functions  for  and  K as  developed  in 

o 

Section  VI.  When  + K2  = 4,  we  deduce  the  approximate  optimum 
Mach  number  M = 0,851  with  the  resulting  value  for  b = 5.799053. 
Hence,  in  both  cases,  the  reduction  in  the  Mach  number  is  relatively 
small  as  compared  to  the  reduction  in  the  weight. 

Singular  Aerodynamics  and  Thrust  Magnitude  Controls  For  Maximum 
Endurance 

We  now  derive  the  complete  solution  for  the  maximum  endurance 
problem.  Again  we  restrict  our  consideration  for  flight  in  the 
stratosphere.  For  flight  in  the  troposphere,  the  approach  is  the 
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same  but  the  derivation  is  more  tedious  and  the  resulting  equations 
are  more  lengthy. 

For  the  problem  of  maximum  endurance,  with  free  range,  m 0, 
and  the  relations  (179)  and  (180)  become 

W ' -Vs  1CD  I2  + “‘c  + <CD  kc)M1  - CL  <209> 

o o M 

"KCLtMkc  " (Kkc)M] } 

and 

Vvv  = -^tccd  - KCJ  + “I  cl)  (210) 

o M 

Rigorously  speaking,  by  eliminating  np^,  Mp^  among  these  equa- 
tions and  the  combined  relation  from  Eqs.  (195)  and  (196) 

Mp  = Mk  (np  + C,)  (211) 

M c n J 

we  have  a relation  for  evaluating  the  optimum  lift  coefficient 

in  terms  of  the  variables  M,  u>,  y and  the  constant  ratio  C^/C^. 

Also  by  replacing  C,  by  p /2KMp_,  in  Eq.  (209)  and  then  taking  its 
L Y M 

derivative,  we  will  have  an  equation  in  which  the  variable  thrust 
magnitude  ht  appears  linearly,  and  thus  can  be  solved  in  terms  of 
the  same  variables  and  constant  ratio.  But  in  doing  it  this  way,  the 
expressions  for  the  aerodynamics  and  thrust  controls  are  unnecessar- 
ily long  and  the  discussion  of  the  way  to  evaluate  the  constant  C^/C^ 
involved  is  not  a simple  matter.  The  simplest  way,  without  compro- 
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mising  the  accuracy  is  to  neglect  the  term  KgnPriY«  which  has  been 
shown  above  to  be  very  small,  in  Eq.  (209).  Then  we  immediately 
have  the  equation  for  the  optimum  lift  control 


**1CL  - ^ CL  + S r2  • 0 

M o 


where  this  time,  we  define  the  function  and  K2  as 


h ' (Kkc>M  - “c 

K2  - 2 + (CD  kc>  + Mkc 


Solving  the  quadratic  equation  (212)  for  C 
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'L  KlM2 


[1 


/ CDK1K2  , 

±Wi  --s-m4 


Km 
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(213) 


(214) 


Maximum  endurance  is  generally  in  the  range  of  Mach  number  where 
compressibility  effect  is  negligible.  Hence  while  K2  is  always  pos- 
itive, is  negative  and  small.  The  two  roots  of  Eq.  (214)  are 
real.  Regarding  the  sign  (+)  in  the  solution,  it  is  seen  that  the 
(-)  sign  must  be  used.  When  is  very  small,  we  have 

CD  K1K2 

C = [1  - (1 2-t—  M4  + . . . ) ] 

KjM  2K<o 

or  approximately  M2 

M k2 


C 


L 


o_ 

2u)K 


(215) 


2 

Since  by  the  equilibrium  cruise  condition  u/M  5 and  at  low  Mach 
number  : 2 we  see  that 
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(216) 


Hence  the  optimum  flight  for  maximum  endurance  Is  conducted  at  near 
maximum  lift-to-drag  ratio.  For  better  accuracy,  the  Eq.  (214)  with 
the  (-)  sign  must  be  used. 

To  have  the  optimum  singular  thrust  control,  we  use  the  optimum 
relation  (127)  to  rewrite  the  equation  (212)  as 

HKK-p  2 - 4K2nwp  pM  + 4K2K2Cd  M3pM2  = 0 (217) 

T ' o 

By  taking  the  derivative  of  this  equation,  using  the  Eqs,  (191) 
with  C2  “ 0 we  have  the  equation  for  the  thrust  magnitude  control  in 
sustaining  flight. 

As  with  the  case  of  maximum  range,  we  first  have 


IKK  p (1  + K^  + K ) 


8K  Wk  . . ,,2 2.  2j 


+ 4K  K2CD  M PM  (3  + 2KM  + 


CD  + K2  )]  - M2(Cd  + KC2)]  - 4K2(KBmMy  - k.nx)p  P,, 

r>  w O 


°M 


(218) 


-2KbKM(MKiPy  - 2KwpM)nPn  + 4K2(2K2CD  M3pM  - up^) 


Mp, 


x{_irtcD  <2  + s > + - T cl]  + nT}  “ 0 


2 2Ku) 


VM 


M 


M 


2 2 


By  using  relation  (127)  and  simplifying  by  4K  p^r  we  have 
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[KKjJT  (1  + Kn  + K^C^  - 4K^Co)Cl  + K2Cd  M^(3  + 21^  + Cd  + K^)  ] 

M 

x^[m  - M2(Cd  + KC2)]  - ZKMC^KgiuMY  - k^x)  + K^w  - KjM^x 


np 


n L 2 


+ ^K2CD.M  - wKCL)  X[M  CD„  (2  + CdJ  + V*1  CL  - 2KwCL  + 


M 


k nx]  = o 

Si 


M 


(219) 


Concerning  the  ratio  np^/Mp^  it  is  obtained  by  solving  the 
Eqs.  (210)  and  (211).  We  have 


np. 


" Mk  (Cn  - KC2  + CT  ) ] (220) 


MpM  UMkc(Cl/C3M  - KgY ) lC3M  c ' Dq  L m2  L' 

As  is  with  the  case  of  maximum  range,  we  shall  neglect  this 
term  and  also  the  term  I^y  in  Eq.  (219).  Then,  solving  for  nx  which 
appears  linearly  in  the  equation,  we  have 


P 

nT  = Q 


(221) 


where  this  time,  we  have 

P = 2mM2KC  (1  - K3  + Kx  ) (CD  + KC2)  + M4(K2  - ^ ) (Cp  + KcJ) 

Mo  o M M o 

+ K2Cd  M4 (K3  + K2)[(Cd  + KC2)  - 2(Mkc/K2)KC2]  (222) 
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and 


kk1m2c^(i  - kx  + ) + k2cd  m2(i  + k2  + k2  ) 

O 


(223) 


'l  / *'2^  »*  \ j.  i ’ *'2  . 

M o M 

As  with  the  case  of  maximum  range,  we  use  this  optimum  thrust  law  in 
the  equation  in  M in  system  (122)  with  n=l  to  have 


2k  C (K.  + K,)M5KC2 
dM  C o 1 * L 

d0  il)Q 


(224) 


This  equation  is  the  exact  equation  for  H.  Again  it  is  seen 
that  the  Mach  number  decreases  slowly.  Combining  this  equation  with 
the  equation  for  w,  we  have 


dw 

dM 


wP 

2Cq  0^  + K2)M5KCl2 
o 


(225) 


If  in  the  expression  (222)  for  P,  we  neglect  the  small  term  Mkc  as 
compared  to  the  drag  coefficient,  and  then  use  the  approximate  sol- 
ution (186),  we  have  the  approximate  equation 


dw  _ bdM 
w M 


(226) 


where  now  b is  defined  as 


2(1  - Kj.  + K ) + K2(K2  - K.)  + K2(Kx  + K2) 

b (k-/k2)  — <227> 


Optimum  Mach  number  for  maximum  endurance  is  low  enough  such  that 
the  effect  of  compressibility  can  be  neglected.  Hence  we  have 
approximately 
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(228) 


b = 3 + 2Mk 

c 


and  we  write  Eq.  (226)  as 


dw 

w 


(3  + 2Mk  )d (Mk  ) 
c c 

(Mkc) 


(229) 


Upon  integrating  from  the  initial  point  to  the  final  point  we  have 

the  approximate  law  for  the  reduction  of  the  Mach  number 

W M 

Log  ^ = 2kc(Mi  - Mf)  + 3 Log  ^ (230) 

Finally,  to  show  the  effect  of  the  altitude,  we  write  the  equa- 
tion for  the  mass  flow 


dw  _ ^c 

de  = " n Q 


(231) 


Now,  if  in  the  expressions  (222)  and  (223)  for  P and  Q,  we  retain 

the  effect  of  Mk  but  using  the  condition  of  equilibrium  cruise  at 
c 

maximum  lift  to  drag  ratio  and  neglecting  the  effect  of  the  Mach  num- 
ber, we  have  after  some  manipulation 


wk  (3  + Mk  ) 
dw  _ c c 

d6  = ~ E (3  + 2Mk  ) 
max  c 


(232) 


Next,  since  M = n1/2w1/,2(K/CD  )1/4,  we  have  the  equation 


o 

1/21/2 


dw 

dd 


wk  [3  + k n ' w 
c c 


(K/cD  )1/4] 

O 


Emax13  + 2kcn1/2w1/2(K/CD  )1/4] 


(233) 
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To  minimize  the  mass  flow  we  minimize  the  right  hand  side  of  this 

equation.  At  low  subsonic  speed  where  is  constant,  this  right 

1/2 

hand  side  is  a homographic  equation  in  n with  negative  derivative. 

1/2 

Hence  it  decreases  as  n increases.  Since  n ■ 2 W^/kpS,  we  should 
use  the  minimum  ambient  pressure,  hence  conduct  the  flight  at  the 
highest  altitude  provided  that  M M^.  Needless  to  say  that  Eq. 

(232)  shows  that  the  effect  of  the  altitude  is  small  and  can  only 
be  detected  by  retaining  the  small  term  Mkc,  that  is  the  second 
order  term.  Also  the  equation  shows  explicitly  that  high  maximum 
lift-to-drag  ratio  favors  saving  in  fuel  consumption. 

Finally,  in  Eq.  (233),  holding  n constant  and  using  constant 
Mach  number,  we  can  integrate  the  equation  to  have  an  approximate 
expression  for  the  time  of  flight 


2E 


V[3  + k n1/2(K/Cn  )1/Aw.1/2) 


max 


Log< 


Wf[3  + kcn1/2(K/CD  )1/4wf1/2] 

O 


(234) 


Again,  in  this  explicit  expression,  it  is  seen  that  the  right 
hand  side  increases  with  n.  Hence,  high  altitude  cruise  favors  en- 
durance. But,  if  the  cruise  altitude  is  taken  in  the  range  where 
the  optimum  Mach  number  is  greater  than  the  effect  of  compressible 

fluid  decreases  E and  the  endurance  decreases  as  the  altitude 
max 

increases. 

A summary  of  the  pertinent  formulas  derived  in  this  section 
and  the  procedure  for  computing  the  optimum  cruise  trajectory  are 
presented  in  Appendix  B, 
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SECTION  VIII 


CONCLUSIONS 

This  report  has  presented  the  results  of  an  extensive  inves- 
tigation of  the  cruise  performance  problem. 

A typical  mission  profile  consisted  of  a climb  at  high  thrust 
setting  to  a certain  optimum  altitude,  followed  by  a cruise  at  var- 
iable thrust  until  a final  prescribed  weight  was  nearly  reached. 

The  flight  terminated  by  a descent  at  reduced  thrust  setting.  For 
long  range  aircraft  the  main  portion  of  the  trajectory  was  the  var- 
iable-thrust portion  called  the  sustaining  arc  along  which,  using 
optimum  control  laws,  significant  savings  in  fuel  consumption  could 
be  made. 

It  was  proposed  to  find  the  variable  thrust  magnitude  and  the 
corresponding  lift  coefficient  in  sustaining  flight  such  that  for  a 
prescribed  final  weight,  either  the  range  or  the  endurance  was  max- 
imized. 

The  problems  were  first  solved  for  the  case  of  constant  alti- 
tude flight  for  both  the  maximum  range  and  the  maximum  endurance 
problems.  It  was  shown  that  the  influence  of  the  Mach  number  on 
aerodynamics  and  engine  characteristics  was  of  prime  consideration 
in  evaluating  the  optimum  Mach  number. 

Next,  to  improve  the  performance,  the  altitude  was  allowed  to 
vary  leading  to  the  concept  of  cruise  climb.  Using  the  approxi- 
mation of  equilibrium  cruise  and  quasi-unaccelerated  flight,  it  was 
shown  as  a first  order  approximation  that  for  both  problems  the  op- 
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timum  Mach  number  was  nearly  constant  and  the  flight  was  at  maximum 
lift-to-drag  ratio.  In  the  maximum  range  problem,  there  existed  an 
optimum  cruise  altitude. 

In  the  maximum  endurance  problem,  the  effect  of  the  altitude  on 
endurance  of  a high  speed  jet-propelled  aircraft  was  small  and  could 
not  be  detected  in  the  first  order  approximation.  The  problems  were 
then  solved  using  singular  optimum  control  theory.  For  both  the 
maximum  range  and  the  maximum  endurance  problems,  the  optimum  lift 
coefficient  and  the  optimum  variable  thrust  magnitude  nr  were 
expressed  explicitly  in  terms  of  the  Mach  number  M and  the  dimen- 
sionless wing  loading  u = 2W/kpS.  Hence,  upon  substitution  into  the 
equations  of  motion,  optimum  trajectory  leading  to  maximum  range  or 
maximum  endurance  could  be  generated  by  numerical  integration.  It 
was  shown  that  in  both  cases  the  optimum  Mach  number  slowly  decreases 
along  the  sustaining  arc.  Explicit  fonnulas  for  evaluating  the  re- 
duction in  the  Mach  number  were  presented.  In  the  case  of  maximum 
range  this  reduction  depended  upon  the  ratio  of  the  final  weight  to 
the  initial  weight  and  the  aerodynamic  characteristics  of  the  air- 
craft. In  the  case  of  maximum  endurance  the  reduction  in  the  Mach 
number  depended  upon  the  weight  ratio  only. 

The  influence  of  the  altitude  on  the  endurance  was  displayed 
explicitly  by  the  rigourous  analysis.  It  was  shown  that  high  max- 
imum lift-to-drag  ratio  at  low  altitude  cruise  favored  endurance. 

A set  of  dimensionless  variables  were  introduced  to  write  the 
equations  of  motion  in  dimensionless  form.  Consequently,  the  re- 
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suits  obtained  were  general  and  apply  to  any  high  speed,  conventional, 
jet  propelled  aircraft. 


In  numerical  applications,  a modelling  of  the  aerodynamics  and 
engine  characteristics  was  necessary.  While  in  the  applications, 
the  specific  fuel  consumption  coefficient  was  assumed  independent 
of  the  Mach  number,  this  effect  was  retained  in  the  pertinent  equa- 
tions derived. 

Two  typical  military  aircraft  were  used  in  the  applications, 
namely  the  transport-tanker  KC-135A  and  the  fighter  F-4C.  Func- 
tional modelling  of  the  aerodynamic  characteristics  of  these  air- 
crafts was  proposed  and  provides  excellent  agreement  with  flight 
test  data  in  the  range  of  Mach  numbers  of  interest. 


APPENDIX  A 


1 

AERODYNAMIC  CHARACTERISTICS 

KC  - 135A 

The  Boeing  KC  - 135A  is  a 4-engine  jet-powered  tanker-trans- 
port and  is  a military  derivation  of  the  Boeing  prototype  model 
367-80. 

The  maximum  lift  coefficient  C^  and  the  buffeting  lift  co- 
max 

efficient  C as  function  of  the  Mach  number  are  given  in  Table  1. 

D 


Table  1 


M 

CL 

max 

CL  ^ 
max 

V* 

0.1 

1.33 

0.013 

1.03 

0.010 

0.2 

1.33 

0.053 

1.03 

0.041 

0.3 

1.31 

0.118 

1.00 

0.090 

0.4 

1.22 

0.195 

0.95 

0.152 

0.5 

1.08 

0.207 

0.88 

0.220 

0.6 

0.99 

0.356 

0.84 

0.302 

0.7 

0.96 

0.470 

0.83 

0.407 

0.8 

0.90 

0.576 

0.66 

0.422 

0.9 

0.68 

0.551 

0.17 

0.138 
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Drag  polars  for  the  KC-135A,  as  function  of  the  Mach  number, 
have  been  obtained  from  a correlation  of  wind  tunnel  tests  and  flight 


test  results.  From  the  drag  polars,  the  values  of  Cp  and  K for 

o 

M J_0.7  are  deduced.  The  experimental  data  and  the  numerical  values 
from  mathematical  modelling,  Eqs.  (A-l)  and  (A-2)  are  presented  in 
Table  2. 


Table  2 


M 

Experimental  Data 

Numerical  Modelling 

0.70 

0.01230 

0.05864 

0.012300 

0.05864 

0.72 

0.012317 

0.05940 

0.74 

0.012338 

0.06056 

0.76 

0.012365 

0.06237 

0.78 

0.01240 

0.06523 

0.012400 

0.06523 

0.80 

0.01245 

0.06991 

0.012450 

0.06992 

0.82 

0.01254 

0.07683 

0.012530 

0.07787 

0.84 

0.01290 

0.09222 

0.012641 

0.09214 

0.86 

0.01350 

0.11958 

0.013500 

0.11961 

0.0123  + 


2.29xl0~*(M  - 0.7) 


(0.8665  - M) 


0.69339 


(A-l 


K 


0.05864 


+ 


7.3xlO~^(M  - 
(0.9742  - M) 


0.7) 

2.884012 


(A-2! 


i 


F - 4C 


The  maximum  lift  coefficient  C and  the  buffeting  lift  co- 
max 

efficient  C as  function  of  the  Mach  number  for  the  fighter  F-4C 
LB 


are  given  in  Table  3 


Table  3 
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For  the  F-4C,  experimental  values  of  CQ  and  K for  M ^ 0. 7 and 

o 

their  corresponding  values  from  mathematical  modelling,  Eqs.  (A-3) 
and  (A-4)  are  presented  in  Table  4 


Table  4 


M 

Experimen 

CD 

o 

tal  Data 

K 

Numerical 

CD 

o 

Modelling 

K 

0.70 

0.01540 

0.155203 

0.015400 

0.155203 

0.72 

0.015516 

0.157578 

0.74 

0.015646 

0.160121 

0.75 

0.01572 

0.161489 

0.015717 

0.161467 

0.76 

0.015792 

0.162867 

0.78 

0.015958 

0.165855 

0.80 

0.01615 

0.169144 

0.016150 

0.169144 

0.85 

0.01675 

0.178884 

0.016806 

0.179288 

0.90 

0.01800 

0.194771 

0.018000 

0.194771 

0.95 

0.02150 

0.229613 

0,021490 

0.229607 

-3, 


0.0154  + —2 77xlO_  (M  ~ Q_‘ 7 ) 


(0.9848  - M) 


0.706112 


(A-3) 


K 


0.155203 


+ 


6.573x10 

(0.983 


2(M  - 0.7) 
m)0.442726 


(A-4) 
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APPENDIX  B 


COMPUTATION  OF  THE  OPTIMUM  TRAJECTORY  IN  SUSTAINED  FLIGHT 


Maximum  Range 


Equation  for  optimum  Mach  number 

«CD  *c\  - 2 

O M 

Definition  of  logarithmic  derivative 

d(Log  X)  _ M dX 
*M  d(Log  M)  X dM 


(B-l) 


(B-2) 


Equation  for  optimum  cruise  altitude 


(B-3) 


Optimum  lift  coefficient 


(B-4) 


First,  compute  the  optimum  Mach  number  from  Eq,  (B-l)  and  then 
the  corresponding  w from  Eq.  (B-3).  It  is  assumed  that  a suboptimum 
law  for  the  climb  portion  of  the  trajectory  has  been  used.  Then 
the  optimum  cruise  altitude  is  reached  when  the  following  condition 
is  satisfied 


(B-5) 
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The  resulting  weight  W = W±  is  the  initial  weight  for  cruise  and  the 
resulting  ambient  pressure  p **  pA  is  the  initial  pressure,  and  hence 
the  initial  altitude  for  cruise. 

The  initial  value  of  n is 


n 


i 


(B-6) 


By  integrating  the  equation  of  motion  from  W = Y = Yi  = 0,  n = 

ni  using  constant  optimum  Mach  number  and  the  lift  control  law  (B-4) 
until  W ” we  obtain  the  maximum  range  in  cruise. 

Constant  Mach  number  at  maximum  lift-to-drag  ratio  implies  that 
the  thrust  law  is 


nx  = 2Cd  M2  (B-7) 

o 

To  improve  performance  in  the  maximum  range  problem  we  shall 
use  the  following  singular  thrust  control 


where 


4(2  - K,  + K,  ) + K2(Kx  + K,  + K0  - K,  ) 


nx  = <2Cd  M ) 
o 


Hi 


2m  y 


K (2  - K + K ) + K,(2  + K,  + K,  ) 

1 XM  z Z Hi 


(B-8) 


K - 1 + (Kk  )v  - Mk 

1 c M c 

K - 1 + (C_  k )M  + Mk 

2 D c M c 

o 


(B-9) 


If  this  law  is  used,  the  optimum  Mach  number  slowly  decreases  along 
the  sustained  arc.  In  this  case,  the  integration  starts  with  an  in- 
itial Mach  number  M^  slightly  larger  than  the  average  optimum  Mach 
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number  computed  from  Eq.  (B-l) . This  value  is  a parameter  to  be 
optimized  for  best  range. 

Maximum  Endurance 


Optimum  lift  coefficient 


(B-10) 


The  altitude  only  has  a secondary  effect  on  endurance.  It  can  be 
taken  arbitrarily  at  low  altitude  where  M <_  Mj  with  M^  being  the 
Mach  number  beyond  which  the  effect  of  compressible  fluid  being 
felt. 

For  the  cruise  arc,  the  integration  starts  with  the  initial 
condition 


w = Wif  n = n±  = 2W1/kpiS,  y * yi  = 0,  M = Mt  (B-ll) 

The  average  value  of  the  optimum  Mach  number  at  any  given  altitude 
is  obtained  from 


(B-12) 


The  Initial  value  M^  of  the  Mach  number  in  Eq.  (B-ll)  is  selected 
slightly  higher  than  the  average  value  of  M obtained  from  Eq.  (B-12), 
and  is  a parameter  to  be  optimized.  It  is  to  be  selected  to  obtain 
the  resulting  best  endurance. 


T 
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The  singular  thrust  control  is  given  by  Eqs.  (221)-(223)  which, 

2 

upon  simplification  by  using  * vC^  /K,  u>  = M and  neglecting 

o 

the  small  term  Mk  , become 
c 


nt 


2(1  - K.  + K.  ) + K,(K  + K,  + K - K,  ) 

(2C  h2  M I JL 

" ; K,  (1  - Kx  + Kx  ) + K2(l  + K2  + K2  ) 


M 


M 


(B-13) 


where 


K1  ' «Vm  - "“c 


K2  ■ 2 + (CD  kc)M  + Mkc 
o 


(B-14) 
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